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ABSTRACT 


This  computer  simulation  investigates  the  ways  in  which  an  embedding 
lattice  influences  a  collision  event.  Low  energy  events  were  studied  with 
a  Born-Mayer  potential  function;  high  energy  with  a  Thomas-Fermi-Firsov 
potential.  Conclusions  were:  (1)  the  lattice  will  increase  ranges  for 
particle  energies  above  400  ev;  (2)  the  apparent  mass  concept  is  not  a 
valid  description  of  events  in  a  lattice;  (3)  lattice  effects  will  signi¬ 
ficantly  modify  the  low  energy  portion  of  the  target  atom  energy  distri¬ 
bution  function;  (4)  there  is  no  evidence  that  a  copper  atom  will  "’rebound88 
from  a  copper  lattice. 

The  writer  wishes  to  express  his  appreciation  for  the  assistance  and 
encouragement  given  him  by  Associate  Professor  Don  E.  Harrison*,  Jr.  of  the 
U.  S.  Naval  Postgraduate  School  in  this  investigation. 
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1.  History. 

Radiation  damage  is  the  overall  term  applied  to  erosion,  disruption*  or 
rearrangement  of  a  crystal  lattice  by  bombardment  of  ions  or  atoms.  Sput¬ 
tering,  a  phenomena  sometimes  observed  in  radiation  damage,  applies  to  the 

process  of  lattice  atoms  being  ejected  or  "knocked  out**  of  the  crystal.  The 

1 

sputtering  process  was  first  observed  in  1852  by  Grove  in  the  cathodes  of 
electrical  discharge  tubes. 

Since  Grove's  observation,  experimental  and  theoretical  physicists  have 
been  attempting  to  develop  a  sputtering  theory  that  will  predict  measurable 
experimental  quantities  as  a  function  of  mass  ratio,  position  in  the  periodic 
table,  energy  of  bombarding  particle  and  angle  of  incidence.  To  date*  no 
such  complete  theory  is  available. 

A  tremendous  volume  of  material  has  been  written  on  the  subject  of  spot- 

2 

tering  since  its  discovery.  In  1955,  Irvine  listed  306  papers  concerned 

with  the  sputtering  process.  These  early  papers  disagree  violently  because 

important  experimental  conditions  were  not  obtained  or  simply  ignored; 

consequently  many  of  the  conclusions  reached  are  contradictory. 

During  this  period  two  important  theories  were  proposed.  Kingdon  and 
3 

Langmuir  ,  in  1923,  developed  a  momentum  transfer  theory  which  agreed  quite 
well  with  experimental  data  they  had  obtained  for  thoriated  tungsten  fila¬ 
ments.  The  "evaporation  theory"  was  proposed  in  1928  by  von  Hippel  and 
4  5 

Blechschmidt  and  improved  by  Townes  in  1944.  It  assumes  that  the  incident 

*W.  R.  Grove,  Phil.  Trans.  Roy.  Soc.  London,  142,  87  (1842). 

2 

M.  M.  Irvine,  Dissertation,  Lehigh  University  (1955). 

3 

K.  H.  Kingdon  and  I.  Langmuir,  Phys.  Rev.  _22,  148  (1923). 

4 

A.  von  Hippel  and  E.  Blechschmidt,  Ann.  Physik  336,  1006  (1928). 

"*C.  H.  Townes,  Phys.  Rev.  65,  319  (1944). 
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particles  heat  a  small  volume  of  the  crystal  to  a  very  high  temperature 

6 

and  subsequently  surface  atoms  evaporate.  More  recently  (1953) 9  Harrison  , 
following  a  suggestion  made  by  Keywell^,  presented  a  statistical  theory  of 
sputtering  which  could  be  placed  somewhere  between  the  two  models  mentioned 
above.  We  will  not  present  details  of  these  theories.  We  must  mention  how¬ 
ever,  that  they  do  not  provide  an  adequate  description  of  the  sputtering 
process. 

The  trend  in  sputtering  theories  has  recently  been  toward  the  considers- 

& 

tion  of  the  individual  collision  processes  within  the  lattice.  Hensehke 

9 

has  used  this  approach  in  his  theory  of  sputtering  and  Harrison  and  Magnuson 
have  also  applied  it,  with  modifications,  to  the  theoretical  study  of  sput¬ 
tering  threshold  energies. 

The  theoretical  work  performed  by  Hensehke  uses  the  momentum  transfer 

3 

concept  initiated  by  Kingdon  and  Langmuir  to  account  for  all  the  experi¬ 
mental  phenomena  observed  at  low  incident  particle  energy.  His  basic  treat¬ 
ment  of  collisions  between  atoms  in  a  lattice  is  based  on  the  assumption 

...that  the  collisions  described  can  be  treated  with  the  general 
principles  of  classical  mechanics,  using  impulsive  forces,  in  a 
manner  similar  to  the  well-known  collisions  with  restitution.^ 


However,  he  did  not  in  all  cases  use  the  masses  of  the  individual  particles 
involved  in  the  equations  pertaining  to  the  two  body  collisions.  Instead 
he  postulated  an  "ef fective"  mass.  If  the  collision  is  between  a  surface 

^D.  E.  Harrison,  Jr.,  Phys.  Rev.  102a  1473  (1956). 

^F.  Keywell,  Phys.  Rev.  87,  160  (1952). 

•  B.  Hensehke,  Phys.  Rev.  106,  737  (1957). 
q 

D.  E.  Harrison  and  G.  D.  Magnuson,  Phys.  Rev.  122s  1421  (1961) 

^See  reference  8,  p.  738. 
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atom  and  a  moving  ion  or  atom  and  is  such  that  the  surface  atom  is  struck 
on  its  "inside"  hemisphere,  according  to  Henschke,  an  effective  mass  is  not 
required  and  the  two  particles  involved  can  be  considered  to  have  their  actual 
masses.  Should  the  collision  of  an  atom  or  ion  be  directed  inward  from  the 
target  lattice  surface,  Henschke  states, 

The  bulk  of  the  target  is  behind  the  struck  atom  and  produces  a 
very  large  ’effective1  mass  compared  to  the  mass  of  the  ion  or 
to  the  mass  of  the  target  atom.H 


The  reason  Henschke  did  not  assume  pure  elastic  collisons  was  stated  as 


Energy  losses  are  due  to  the  fact  that  the  target  atom  is  coupled 
rather  strongly  with  the  atoms  of  the  lattice.  Before  the  moment 
of  highest  compression  is  reached,  the  ion  and  the  struck  atom 
exchange  energy  with  the  neighboring  atoms  of  the  lattice.  Debye 
waves  are  thus  excited  and  dissipated  irreversibly  into  the  lat¬ 
tice. 


He  goes  on  to  state. 


...the  final  step  in  each  sputtering  process  at  any  angle  of  inci¬ 
dence  of  the  ion  can  be  generally  described  as  a  collision  of  the 
ion  with  an  upper  surface  atom,  in  which  this  atom  is  hit  on  its 
inside  hemisphere  so  as  to  obtain  an  impulse  with  a  component  in 
the  direction  of  the  outward  normal  to  the  surface.  If  the  energy 
transferred  in  this  direction  to  this  target  atom  by  the  impart 
of  the  ion  is  equal  to  or  greater  than  the  heat  of  vaporization, 
with  which  the  atom  is  assumed  to  be  bound  to  the  crystal  lattice 
plane,  then  this  atom  is  ejected  in  the  collision.^ 


Early  study  in  radiation  damage  was  precipited  by  the  advent  of  the 
nuclear  reactor.  Recently  however,  with  the  imminent  possibility  of  thermo¬ 
nuclear  power  production  and  ion  propulsion  engines  for  space  vehicles,  other 
practical  applications  of  the  processes  involved  in  radiation  damage  have 


^E.  B.  Henschke,  Phys.  Rev.  121,  1290  (1961). 
12 

See  reference  8,  p.  738. 

13 

See  reference  8,  p.  737. 
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become  extremely  important, 

14 

Robinson  and  his  co-workers  have  made  theoretical  studies  of  the 
ranges  in  solids  of  atoms  having  energies  from  1  to  100  kev  using  digital 
computer  techniques  on  the  basic  assumption  that  the  moving  atom  loses  its 
energy  through  repeated  binary  elastic  collisions  with  atoms  of  the  solid. 

The  masses  used  in  these  calculations  are  the  true  masses  of  the  interact¬ 
ing  particles,  but  they  note  that  their  assumption  is  certainly  not  valid 
below  100  ev,  Gibson  et.  al,^  have  also  taken  advantage  of  the  speed  avail¬ 
able  with  modern  digital  computers  in  the  study  of  radiation  damage.  They 
do  not  assume  binary  collisions  but  instead  use  iteration  techniques  employ¬ 
ing  Newton* s  equation  of  motion  to  solve  the  complex  many  body  problem.  The 
initial  success  of  these  programs  have  encouraged  their  originators  to  ex¬ 
plore  further  the  possibilities  of  these  techniques.  Although  the  computer 
programs  are  designed  to  study  radiation  damage,  the  basic  principles  involv¬ 
ed  are  also  important  in  sputtering. 

The  work  just  discussed  is  very  sensitive  to  the  mathematical  form  of 
the  interatomic  potential.  The  Born-Mayer,  Bohr  (screened  Coulomb),  and 
the  Thomas-Fermi-Firsov  are  some  of  the  potentials  used  in  these  calcula¬ 
tions.  The  ability  of  any  of  the  potential  functions  to  describe  a  physical 
situation  is  highly  dependent  upon  the  energy  range  under  consideration  and 
the  specific  atoms  or  ions  involved.  No  single  potential  has  yet  been  de¬ 
vised  that  satisfactorily  represents  the  interaction  under  all  circumstances. 

14 

0.  S.  Oen,  D.  K.  Holmes,  and  M.  T.  Robinson,  Jour  Appl.  Phys.  34, 

302  (1963). 

^J.  B.  Gibson,  A.  N.  Goland,  M.  Milgram,  and  G.  H.  Vineyard,  Phys. 

Rev.  120,  1229  (1960) 
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An  obvious  unresolved  question  is  inherent  in  these  theories 0  The 
many  body  approach  and  the  two  body  assumption,  with  a  minimum  energy  limita¬ 
tion,  appear  to  be  reasonable  approaches  to  the  problem0  Both  methods  des¬ 
cribe  possible  events  in  a  lattice,  but  are  they  compatable?  The  general 
use  of  the  two  body  approach  is  desirable  because  of  its  simpliclty0  The 

problem  to  be  studied  is:  when  can  the  two  body  collision  be  assumed,  if 

8 

at  all,  and  should  an  "apparent”  mass,  as  proposed  by  Henschke  ,  be  employed? 


r 
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2,  Object ive. 

We  intend  to  investigate  theoretically  the  energy  transfer  process 

between  an  incoming  copper  atom  and  a  copper  lattice.  The  many  body  ap~ 

15 

proach,  using  computer  techniques  similar  to  Gibson  efc.  aL,  will  be 
utilized  and  then  compared  with  the  simple  two  body  solution.  We  hope 
the  results  will  provide  some  quantitative  answers  to  the  questions  posed 
in  the  preceding  section. 
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3,  Procedure 


The  representative  crystal  lattice  consists  of  a  cube  of  63  atoms 
arranged  in  a  face  centered  cubic  structure  *  This  arrangement  gives  a 
total  of  63  atoms  *  These  atoms  are  free  to  move*,  if  disturbed.  In  an 
effort  to  simulate  a  larger  lattice,  the  movable  63  atoms  are  surrounded 
on  all  sides  by  stationary  immovable  atoms  positioned  as  a  continuation  of 
the  crystal.  This  procedure  also  gives  the  bombarded  face  some  character- 
istics  similar  to  a  binding  energy. 

In  order  to  simplify  the  calculat ions *,  we  defined  a  quantity  called 
the  "lattice  unit".  One  lattice  unit  is  equal  to  half  the  length  of  a 
cube  side.  For  copper,  the  length  of  a  cube  side  is  3,614A&  therefore 
one  lattice  unit  is  1,807A.  A  three  dimensional  Cartesian  coordinate 
system  is  used  for  reference  with  the  origin  placed  with  a  (100)  face  of 
the  movable  63  atom  "core"  parallel  to  the  x-z  plane  at  y-0,  Any  penetra¬ 
tion  into  the  lattice  core  from  this  plane  is  in  the  4*  y  direction.  The 
immovable  atoms  surrounding  the  core  are  in  the  x-z  plane  at  y--l  and  y-59 
the  x~y  plane  at  z-0  and  z=6,  and  the  y-z  plane  at  x^Q  and  x-6,  For  conven¬ 
ience  all  atoms  are  numbered.  Number  one  is  the  "bullet"*,  numbers  2-64  are 
the  "core"  and  65-172  are  the  immovable  atoms,  (see  Fig,  1) 

The  program  depends  upon  a  potential  function  subroutine;  so  potential 
modification  is  accomplished  without  any  changes  in  the  main  program.  The 
Bom-Mayer  and  the  Thomas-Fermi-Firsov  have  been  used  in  the  present  work 
(see  Appendix  VI  for  exact  forms  and  constants).  All  calculations  are 
with  the  potential  in  "eroded"  form.  To  "erode"  a  potential^,  the  value  of 
the  potential  at  the  nearest  neighbor  separation  is  subtracted  from  the 
value  of  the  potential  at  smaller  distances.  The  eroded  form  makes  the 
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potential  zero  at  distances  greater  than  nearest  neighbors. 


The  interaction  between  atoms  is  represented  by  an  eroded  repulsive 
potential.  Also,  the  entire  lattice  is  initially  at  absolute  zero  since 
no  vibrational  energy  is  simulated.  These  two  approximations  imply  that 
the  lattice  has  no  potential  or  kinetic  energy  before  interaction  with  the 
bullet  and  all  the  energy  in  the  lattice  at  any  time  thereafter  is  derived 
from  the  bombarding  atom. 

We  realize  that  each  atom  in  the  lattice  is  contained  in  some  type  of 
"potential  well",  but  we  are  measuring  the  energy  from  the  ground  state 
level  of  this  "well"  rather  than  from  the  true  zero  of  potential  at  infinite 

separations.  Thus,  before  interactions  with  the  bullet,  our  lattice  is  not 

15 

held  together  by  the  surrounding  stationary  atoms  as  in  the  Gibson  et.  al.“ 
model  because  the  erosion  also  has  removed  all  forces.  Forces  appear  only 
when  atoms  move  from  their  equilibrium  positions. 

The  use  of  an  eroded  potential  does  not  provide  as  good  a  model  as  the 
one  used  by  Gibson  et.  al,,  but  it  reduces  the  computer  time  required  for 
calculations. 

The  "binding  energy"  created  by  the  layer  of  fixed  atoms  covering  the 
front  face  is  probably  much  too  large,  but  it  does  approximately  simulate 
the  behavior  of  the  next  layer  of  atoms. 

IS 

The  calculations  which  move  atoms  are  similar  to  those  of  Gibson  et0  al* 
Newton® s  equation  of  motion  can  be  rearranged  to  give  the  change  in  velocity 


of  a  body  acted  on  by  an  unbalanced  force*  ( 


The  change  in  velocity  can  be  related  to  a  change  in  position  if  an  average 


velocity  is  assumed 


used  in  the  first  equation  is  an  average  force  calculated  by  a  double  itera¬ 
tion  procedure  as  follows:  (1)  assume  an  atom  at  position  1  with  velocity  1 
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(2)  calculate  the  total  force  on  the  atom  as  a  result  of  all  the  other  atoms 
in  the  lattice  (this  means  normally  only  about  8-10  nearest  atoms  because 
the  potential  is  eroded)  (3)  call  this  calculated  force*,  force  18  and  use  the 
equation  of  motion  to  move  the  atom  to  a  temporary  position*,  position  2C 

(4)  now  repeat  the  force  calculations  for  position  2,  call  this  force  20 

(5)  go  back  to  position  1,  and  use  the  average  of  force  1  and  force  2  to 
move  the  atom  to  a  new  position,  position  3,  Procedures  1  through  5  consti- 
tute  one  "time  step!’.  Forces  are  eroded  in  the  same  manner  as  potentials 
and  are  calculated  by  a  subroutine  based  on  the  partial  derivative  of  the 
potential  function  with  respect  to  distance. 

The  basic  value  of  AT  is  dependent  upon  the  original  energy  of  the 
bombarding  atom.  A  T  is  the  time  in  seconds  required  for  the  incoming  atom 
to  traverse  one  lattice  unit.  The  atom  will  lose  energy  as  it  interacts 
with  the  lattice  but  AT  remains  unchanged.  To  add  flexibility  to  the 
program,  the  basic  value  of  AT  can  be  multiplied  by  any  desired  factor 
(time  step  multiplier).  A  T  is  established  at  the  beginning  of  the  program 
and  is  constant  for  the  duration  of  the  calculations. 

The  bombarding  copper  atom  is  originally  located  in  the  plane y  -  ->J~2 
and  given  a  velocity  in  the  +  y  direction  only.  This  requires  the  bombarding 
particle  to  approach  perpendicular  to  the  (100)  face  of  the  movable  lattice 
core.  The  immovable  atoms  in  the  plane  (y  =  -1)  that  cover  the  bombarded 
face  do  not  interact  with  the  bullet. 

An  impact  area  was  chosen  on  the  (100)  face  so  that  all  points  in  this 
area  would  be  representative  of  any  point  in  a  (100)  plane.  An  impact  point 
is  defined  as  the  location  on  the  (100)  face  toward  which  the  bombarding  atom 
is  directed.  It  may  or  may  not  actually  pass  through  this  point.  By  moving 
in  1/10  lattice  unit  increments,  36  impact  points  in  the  impact  "triangle" 
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were  assigned.  (see  Fig.  2  and  3). 

For  each  case  the  two  body  problem  is  solved  after  the  many  body 
problem  with  the  same  double  iteration  scheme.  Lattice  atom  number  eight 
was  assigned  as  the  "target"  in  the  many  body  lattice  problem  so  that  its 
motion  could  be  compared  with  the  simple  two  body  problem.  The  geometrical 
relationship  and  physical  constants  of  the  two  body  problem  correspond  to 
those  of  the  incoming  atom  and  "target"  in  the  crystal. 

The  program  was  written  in  FORTRAN  language  for  use  on  a  CDC  1604 
computer.  The  input  consists  of  the  incoming  particle  mass  and  energy  and 
the  time  step  multiplier.  From  these  inputs  the  bullet  is  assigned  a  loca¬ 
tion  and  an  appropriate  velocity  component  in  the  +  y  direction.  Also  cal¬ 
culated  is  the  actual  value  of  AT. 

The  next  section  of  the  program  assigns  coordinates  to  all  atoms  in 
the  crystal,  movable  and  immovable,  gives  each  atom  a  number  and  sets  their 
velocities  equal  to  zero.  If  the  distance  between  any  two  atoms  is  less  than 
sT l  lattice  units,  the  force  between  these  two  atoms  is  calculated  and  stored 
in  memory.  This  is  done  for  each  atom  and  the  summation  of  all  the  forces 
calculated  gives  the  resultant  force  on  each  particle. 

The  change  in  the  velocities  of  all  movable  atoms  as  a  result  of  un¬ 
balanced  forces  is  calculated;  then  the  original  velocities  of  all  particles 
and  the  calculated  change  in  velocity  are  used  to  move  atoms  to  temporary 
locations.  The  calculations  to  obtain  resultant  forces  are  repeated  for  the 
temporary  locations,  old  and  new  forces  are  averaged  and  used  to  move  all 
atoms  to  n$w  locations.  Forces  and  velocities  are  calculated  as  vector 
components  using  the  coordinate  system  established  for  the  lattice. 

After  each  movement  of  all  the  atoms,  the  kinetic  energy  of  the  tar* 
get  is  calculated  and  compared  to  its  former  value.  The  procedure  is  repeated 
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until  the  kinetic  energy  of  the  target  reaches  a  maximum.  The  potential 
and  kinetic  energy  of  each  movable  atom  is  then  calculated.  The  total 
energy  of  the  entire  lattice  including  bullet  is  also  found  and  compared  to 
the  original  "bullet"  energy.  This  comparison  gives  a  reasonable  measure  of 
the  accuracy  achieved.  The  two  body  problem  is  solved  in  the  next  section 
using  the  same  iteration  technique.  Appendices  I-V  contain  a  detailed  ex¬ 
planation  of  the  program. 

For  each  "run",  the  computer  will  solve  the  many  body  and  two  body 
problem  36  times,  corresponding  to  the  36  different  impact  points.  The  form 
of  the  printed  output  for  each  impact  point  is  shown  in  Fig.  4  and  5.  The 
"triangular"  output  as  seen  in  Fig.  6  indicates  how  two  chosen  parameters 
vary  over  the  impact  area. 

The  computer  running  time  for  each  of  the  36  impact  points  was  2-4 
minutes,  the  time  varying  with  the  value  chosen  for  AT.  We  originally 
assumed  that  the  accuracy  of  the  results  would  increase  for  smaller  and 
smaller  values  of  AT.  We  found  that  the  value  of  AT  which  will  give  the 
minimum  error  in  total  energy  appears  to  be  a  complicated  function  of  the  time 
step  multiplier,  the  incoming  bullet  energy,  and  the  impact  point.  By  trial 
and  error  methods,  we  found  that  below  100  ev  for  the  Born-Mayer  potential, 
a  time  step  multiplier  of  around  0.09  would  produce  reasonable  errors  at  all 
impact  points. 
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4.  Potentials 


The  Born-Mayer  potential  is  an  exponential  function  of  the  internuelear 
separation  with  constants  which  may  be  obtained  from  experimentally  measur¬ 
able  elastic  moduli.  The  constants  used  in  the  program  are  the  same  as 
those  of  Gibson  et.  al.^,  with  special  emphasis  on  their  potential  number 
two.  Near  equilibrium  separations,  the  Born-Mayer  potential  is  thought  to 
be  an  adequate  approximation.  We  expect  the  potential  to  fail  for  energies 
above  a  few  hundred  electron  volts. 

The  Thomas-Fermi-Firsov  (TFF)  potential  is  the  result  of  theoretical 

17 

work  by  Firsov  based  on  the  Thomas-Fermi  model  of  the  atom.  It  is  a 

screened  Coulomb  potential  with  a  more  complicated  screening  function  than 

that  used  in  the  Bohr  potential.  The  TFF  is  satisfactory  for  intermediate 

separations.  This  places  it  between  the  Bohr  potential  used  for  small  separa 

18 

tions  and  the  Born-Mayer  potential.  Abrahamson  and  Hatcher  state  that  the 

Thomas -Fermi  approximation  becomes  unreliable  when  the  intermuelsar  distance 

exceeds  /VIA.  The  form  of  the  TFF  potential  and  the  appropriate  constants 

used  in  the  program  are  the  same  as  those  used  in  computer  programs  at  the 

19 

Oak  Ridge  National  Laboratory 

The  mathematical  forms  of  the  potential  and  force  functions  as  they  ap¬ 
pear  in  the  program  are  derived  in  Appendix  VI. 


^See  reference  15,  p.  1233. 

170.  B.  Firsov,  J.  Exptl.  Theoret.  Phys.  (U.S.S.R.)  32,  1464  (1957). 
translation:  Soviet  Phys.  -JETP  _5,  1192  (1957). 

^A.  a.  Abrahamson  and  R.  D.  Hatcher,  Phys.  Rev.  121 ,  159  (1961). 

19 

R.  T.  Robinson  (private  communication).  The  author  would  like  to 
thank  Dr.  Robinson  for  the  use  of  his  results  prior  to  their  publication. 
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5.  Results 


A.  General  description  of  events  in  the  lattice. 

The  number  of  atoms  disturbed  in  the  lattice  is  a  function  of  two 
variables,  the  number  of  time  steps  the  computer  performs  and  the  impact 
point  of  the  bullet.  In  general,  the  lattice  is  more  disrupted  by  impact 
points  in  the  area  near  x  =  3.0,  z  «  4.0,  and  the  lattice  is  least  af- 
fected  by  head  on  collisions  with  the  target.  No  more  than  one  half  of  the 
atoms  were  ever  disturbed  before  the  target  reached  a  maximum  in  kinetic 
energy. 

The  impact  area  is  not  symmetrically  located  with  respect  to  the  sur¬ 
rounding  immovable  atoms  but  for  appropriate  points  in  the  impact  area, 
excellent  symmetry  of  displacement  was  observed  for  the  movable  atoms.  This 
indicates  that  the  size  of  the  model  is  adequate  for  present  purposes. 

Because  of  the  limitation  placed  on  the  number  of  time  steps,  only  a 
few  atoms  have  any  significant  kinetic  energy  when  the  interaction  is  stop¬ 
ped.  These  are  normally  the  target,  bullet,  and  the  atoms  directly  behind 
the  target  and  bullet  in  the  lattice. 

Although  our  primary  purpose  was  to  compare  lattice  interactions  with 
the  two  body  interaction,  we  altered  the  program  so  that  it  would  run  a 
predetermined  number  of  time  steps  in  order  to  observe  "chains"  and  "chan¬ 
nels".  Impact  points  at  x  =  3.0,  z  =  4.0  and  x  =  3.0,  z  =  3.0  correspond 
to  the  beginning  of  (100)chains.  These  are  not  "close  packed"  chains  but 
they  transfer  energy  reasonably  well,  especially  above  100  ev.  Impact  point 
x  =  2.5,  z  =  3.5  is  the  beginning  of  a  (100)  channel.  The  (100)  channel  is 

■  m  v~ 

not  as  "wide  open"  as  a  (110)  channel  and  a  100  ev  bullet  does  not  travel 
more  than  C*J  3A  into  the  lattice. 
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B.  Comparison  of  lattice  interactions  with  the  two  body  interaction. 

A  strict  comparison  of  interactions  in  the  lattice  to  the  two  body 
problem  is  not  possible,  but  certain  limited  correlations  can  be  made.  The 
difficulty  arises  because  the  end  of  the  interaction  in  the  lattice  cannot 
be  defined  explicitly.  There  are  several  criteria  available  that  may  be 
used  in  an  attempt  to  compare  the  interactions.  If  the  two  body  interaction 
was  allowed  to  continue  until  the  potential  between  atoms  was  zero,  for 
comparison  purposes  the  lattice  interaction  could  procede  until:  1)  the 
potential  between  target  and  bullet  is  equal  to  zero  2)  the  total  energy 
(kinetic  plus  potential)  of  the  target  reaches  a  maximum  or  3)  the  kinetic 
energy  of  the  target  reaches  a  maximum.  All  of  the  criteria  apply  to  the 
same  physical  situation  in  the  two  body  problem  but  are  not  equivalent  in 
the  lattice  interaction. 

We  chose  the  maximum  kinetic  energy  of  the  target  as  the  comparison 
stopping  point  because  this  represents  the  time  when  the  target  and  bullet 
have  almost  ceased  to  interact  and  are  just  beginning  to  interact  with  other 
lattice  atoms. 

A  fourth  criteria  was  suggested  after  the  present  results  had  been 
obtained.  When  the  kinetic  energy  of  the  target  in  the  lattice  reaches  a 
maximum,  find  the  distance  between  the  target  and  bullet  and  allow  the  two 
body  interaction  to  proceed  until  the  distance  between  atoms  is  equal  to  the 
separation  in  the  lattice.  This  criteria  is  possibly  better  than  the  one 
chosen  but  at  the  present  time  it  has  not  been  utilized. 

The  geometrical  relationships  of  the  target,  bullet,  and  lattice  intro¬ 
duce  complications  into  the  comparison  scheme.  The  motion  of  the  target  in 
the  lattice  is  not  restricted  isotropically  (see  Fig.  7).  If  the  inter¬ 
action  with  the  bullet  moves  the  target  toward  atom  number  six  (shaded  area), 
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the  movement  of  the  target  is  restricted  more  than  if  the  target  motion  is 
initially  toward  the  cross-hatched  area.  Restriction  of  the  target  motion 
lengthens  the  interaction  time  between  bullet  and  target,  therefore  the 
energy  transferred  to  the  target  is  greater  than  in  the  two  body  problem® 

This  effect  is  noticable  only  at  impact  parameters  of  0.5  lattice  units  or 
greater.  At  impact  distances  smaller  than  this  the  interaction  occurs  so 
rapidly  that  the  target  does  not  move  a  significant  amount  and  the  effect  is 
reduced*  The  kinetic  energy  transferred  at  distances  smaller  than  0.5  lattice 
units  is  less  than  that  transferred  in  the  two  body  case  because  there  is 
absorption  of  energy  by  the  remainder  of  the  lattice.  The  geometry  of  the 
bombarded  face  is  not  the  only  factor  that  effects  energy  transfer.  The 
atoms  behind  the  target  in  the  lattice  also  play  an  important  role  but  their 
effect  is  not  as  immediately  obvious  as  those  just  discussed. 

Figure  8  illustrates  some  of  the  points  just  mentioned.  Along  the  line 
A  B  (x  =  3.0,  z  =  3. 0-4.0)  the  target  behavior  in  the  lattice  is  similar  to 
the  two  body  interaction.  These  impact  points  cause  the  target  to  move  in 
the  negative  z  direction  toward  a  "hole1'  in  the  lattice  face.  At  point  D  on 
the  line  A  C  ,  the  effect  of  restricted  target  motion  becomes  obvious,  and 
also  in  the  area  enclosed  by  the  dashed  line.  The  geometrical  effects  of 
the  lattice  are  fairly  predictable  and  no  unusual  phenomena  are  observed. 

For  a  static  system  of  two  equal  mass  atoms  with  a  conservative  repul¬ 
sive  force  between  them,  the  total  potential  energy  of  the  system  can  be 
halved  and  the  result  assigned  as  the  energy  of  the  atom.  If  the  two  atoms 
are  allowed  to  move,  the  total  energy  (potential  plus  kinetic)  of  each  will 

»  i 

be  constant  and  equal  to  the  original  potential  energy.  In  a  system  composed 
of  three  or  more  atoms,  this  is  generally  not  true.  For  any  system  where  the 
atoms  are  not  static,  the  process  of  assigning  an  energy  to  any  one  atom  is 
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no  longer  possible  (there  is  one  exception).  This  is  the  situation  in  the 
lattice. 

The  potential  energy  of  an  atom  in  the  lattice  is  defined  in  the  program 
as  half  the  potential  associated  with  its  position.  This  definition  is  merely 
a  convenience. 

The  area  of  the  impact  triangle  from  point  A  to  the  arc  E-F  (see  Fig.  8) 
is  least  effected  by  the  geometry  of  the  lattice  and  closer  inspection  of 
these  interactions  is  justified.  The  kinetic  energy  transferred  to  the 
n-body  target  is  less  than  that  of  the  two  body  target.  In  percentages, 
the  kinetic  energy  transferred  becomes  less  as  the  bullet  energy  is  decreased, 
(see  Figs  8  through  11).  In  each  case  the  target  also  acquires  some  poten¬ 
tial  energy  (see  Figs.  12  through  15).  As  stated  earlier,  it  is  difficult  to 
assign  a  specific  potential  energy  to  the  target.  We  do  know  the  limits 
involved,  i.e.  the  target  could  eventually  receive  none  of  it,  or  twice  the 
indicated  value.  Either  limiting  situation  is  unlikely. 

If  we  make  the  assumption  that  one-half  of  the  potential  "belongs" 
to  the  target,  then  the  total  energy  of  the  target  is  approximately  equal  to 
the  energy  transferred  in  the  two  body  problem  (within  27o)  for  energies  great¬ 
er  than  50  ev. 

Bullet  behavior  for  the  area  under  discussion  is  also  very  similar  to 
the  two  body  problem  (see  Figs.  16  through  19).  For  the  remainder  of  the 
impact  triangle,  the  bullet  is  strongly  effected  by  other  atoms  in  the  lat¬ 
tice. 

Three  runs  were  made  with  the  TFF  potential  at  energies  of  5,  10,  and 
30  kev  (Figs.  20  through  25).  Since  the  TFF  potential  is  unreliable  beyond 
1A,  the  section  ABC  in  Figs.  20  through  25  is  the  only  part  of  the  impact 
area  that  can  be  considered.  This  section  includes  all  impact  points  that 
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result  in  a  closest  point  of  approach  of  1A  or  less.  We  noted  the  same 

general  agreement  with  the  two  body  interaction  as  found  for  the  Born~Mayer 

potential.  The  only  apparent  difference  is  a  shift  in  the  energy  scale. 

For  a  head-on  collision  in  the  lattice  there  is  very  little  geometrical 

effect,  and  an  "effective"  mass  can  be  described.  At  25  ev,  the  mass  of  the 

target  is  apparently  about  2.5  times  the  bullet  mass,  but  the  energy  trans- 

fer  is  still  approximately  80%.  For  a  500  ev  bullet,  the  effective  mass 

increase  is  negligible.  We  used  the  program,  as  modified  to  observe  chains 

and  channels,  and  could  find  no  evidence  of  bullet  recoil  which  could  be 

8 

attributed  to  an  effectively  heavy  target  as  proposed  by  Henschke  .  Although 
an  effective  mass  concept  for  the  head-on  collision  is  possible,  it  does  not 
properly  describe  the  subsequent  motion  of  either  target  or  bullet.  If  the 
bullet  motion  results  in  a  glancing  hit  with  the  target,  no  single  effective 
mass  can  be  assigned  to  the  target  because  it  is  a  function  of  the  impact 
parameter.  An  average  effective  mass  is  not  applicable  because  the  geometri¬ 
cal  effects  of  the  lattice  on  the  direction  and  energy  of  the  recoil  atom  are 
far  more  significant  than  the  mass  of  the  target  in  a  glancing  hit. 

These  results  apply  to  a  collision  anywhere  in  the  lattice  and  are  not 
limited  to  the  interaction  of  surface  atoms  with  incoming  particles.  The 

g 

results  indicate  that  Henschke' s  apparent  mass  concept  and  the  rebound 
phenomena  associated  with  it  are  not  a  good  description  of  collision  events 
in  a  lattice. 


17 


6.  Conclusions 


A,  The  Bullet • 

For  energies  above  40  ev  in  the  impact  area  AEF  (see  Figs  16  and  17)* 
the  n-body  bullet1 s  kinetic  energy  after  interaction  is  essentually  the 
same  as  the  two  body  bullet.  Between  40  and  400  ev  for  impact  points  out¬ 
side  the  area  AEF,  the  n-body  bullet  has  considerably  less  energy  after  inter¬ 
action  than  the  two  body  bullet.  Above  400  ev,  the  n-body  and  two  body 
energy  transfers  agree  within  37o  for  all  points  in  the  impact  area. 

The  angular  behavior  of  a  bullet  in  the  lattice  is  much  more  complicated. 
Even  at  energies  above  400  ev,  the  scattering  angle  in  the  n-body  problem  is 
affected  by  the  geometry  of  the  lattice.  In  general,  the  scattering  angle  in 
the  lattice  is  smaller  than  that  of  the  two  body  problem.  It  is  possible  that 
the  lattice  is  attempting  to  focus  the  energy  into  preferred  directions.  Our 
bullets  in  the  (100)  direction  appear  to  focus  in  the  (100)  direction*  especial¬ 
ly  at  higher  energies. 

The  cummulative  effects  of  energy  transfer  and  lattice  geometry  will 
affect  the  range  of  energetic  atoms.  The  reduced  scattering  angle  in  the 
n-body  model  for  atoms  with  energies  above  400  ev  should  lead  to  ranges  that 
exceed  those  found  by  two  body  approximation  methods.  Below  400  evs  the 
bullet  loses  more  energy  than  in  the  two  body  case;  so  the  effects  may  cancel* 
or  perhaps  ranges  calculated  by  two  body  methods  may  be  too  large, 

B.  The  Target. 

In  the  area  AEF  (see  Figs.  8,  9,  12,13),  the  n-body  energy  transferred 
to  the  target  is  the  same  (within  2%)  as  that  in  the  two  body  problem  for 
energies  above  50  ev.  For  the  remainder  of  the  impact  area,  the  n-body 
target  receives  more  energy  than  the  two  body  target  for  all  energies  above 
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25  ev.  Isolated  points  at  various  energies  may  not  conform  to  this  rule  but 
the  effect  is  always  present  for  impact  points  near  the  apex  of  the  impact 
area  (x=2.5,  z=3.5).  Near  the  apex,  if  the  bullet  energy  is  large  (300  ev 
or  greater),  the  energy  transferred  to  the  target  is  a  very  small  fraction 
of  the  total  for  both  the  n-body  and  two  body  problems.  However,  the  energy 
transferred  to  the  target  in  the  two  cases  can  differ  by  as  much  as  75%  of  the 
transferred  energy. 

As  a  consequence  of  this  phenomena,  the  low  energy  portion  of  the  n- 
body  target  atom  energy  distribution  (assuming  more  than  one  collision  has 
occurred)  will  be  much  higher  than  that  expected  with  a  two  body  collision 
assumption  for  lattice  interactions.  Precise  measurements  of  energy  trans¬ 
fer  for  large  impact  parameters  do  not  appear  to  be  warranted  in  view  of  the 
results  obtained. 

Unfortunately,  the  face  centered  cubic  lattice  structure  was  incorporate 

20 

ed  in  the  program  before  we  learned  of  Veksler's  experimental  work  with 
molybdenum  targets  (bcc),  but  certain  qualitative  comparisons  are  still 
possible.  Veksler  has  interpreted  the  lattice  behavior  in  terms  of  the  ef~ 
fective  mass  concept  but  his  general  conclusions  are  consistent  with  our 
work.  Our  model  gives  good  evidence  to  support  Veksler 9 s *  posit ion  that  the 
pair  collision  model  using  elastic  spheres  is  not  acceptable. 

Two  body  approximations  inherently  imply  that  certain  information  is 
available  from  preliminary  n-body  calculations.  This  work  attempts  to  answer 
certain  questions  about  interactions  and  cross  sections.  As  anticipated^  the 
results  are  not  definitive  but  they  do  indicate  sensitive  areas  which  require 
further  examination. 

20 

V.  I.  Veksler,  J.  Exptl.  Theoret.  Phys.  (U.S.S.R.)  42,  325  (1962). 
translation:  Soviet  Physics  -  JETP  lf>,  222  (1962). 
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. 


small  circles  indicate  fixed  atoms. 
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X=3. 0 

Fig.  2.  FACE  (100)  of  LATTICE  CORE. 

Numbers  indicate  atoms;  shaded  portion  is  impact  area 
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Fig.  4.  TYPICAL  COMPUTER  OUTPUT  for  an  impact  point 


<N 

O 

♦ 

UJ 

Jf 

o 


H 

■si 


O 

O 

O 

o 

o 

• 

o 

o 


CM 

O 

♦ 

UJ 

in 

O' 

• 


o 

Cp 

o 


UJ 

-J  II 
-J  c*. 
3  OJ 
OC 

•J 

O  CL 
z 

<•  h 


g- 

X 


I—  CL 
h-  Ui 
<  ►- 
in 

oc  UJ 
uj  x 
a  *-» 


& 


G 

G 


O 

a 


N 

>- 

O 

OC 


o  d 
oc  o 

UJ  O' 

• 

oc  m 

UJ 

>•  H 
<  NJ 

x  \ 

i 

oc  ©. 
o  o 

■3  o 


II 

X 


< 

u 

o 


o 

z 

►H 

X 

o 

o 

z 


3 

03 


II 


> 

*3 

II 

o 

z 

•— 

1 

UJ 

o 

o 

1 

UI 

* 

r- 

CO 

oc 

o 

m 

z 

H- 

o 

CT' 

UJ 

UJ 

• 

• 

-J 

cc 

-J 

G 

II 

g 

CO 

X 

M 

o 

o 

1 

a 

»— 

o 

K 

< 

o 

r- 

> 

O' 

-L-» 

g 

-JO 

UJ 

o 

(X 

O 

a 

6' 

o 

00 

o 

UJ 

V 

oo 

V 

• 

z 

O- 

o. 

K-  O 

X 

CM 

UJ* 

6 

o 

CJ 

od 

r_ 

oco 

N.  • 

IN 

CT' 

in 

-fe* 

a 

05 

2 

CM 

m 

2 

> 

X 

II 

• 

o 

Im 

V 

*- 

UJ 

UI 

03 

r— 

o 

O 

H 

UJ 

_l 

o 

•— 

CM 

V-i 

o 

o 

a 

03 

o 

• 

oc 

oc 

oc 

o 

o 

v 

• 

o 

H 

z 

< 

-J 

oc 

V 

II 

in 

& 

CM 

UJ 

I- 

UJ 

a 

< 

CO 

m 

X 

> 

-J 

uo 

UJ 

II 

m 

»— 

UJ 

> 

o 

o 

o 

z 

H 

o 

ox 

o 

M 

II 

UI 

r- 

G 

& 

L- 

L- 

OG 

K 

L- 

o 

O' 

oc 

o 

oc 

cox 

•— 

»- 

UI 

• 

O 

a 

O' 

o 

• 

< 

o 

y- 

u 

« 

o 

oc 

ox 

•“ 

o 

-J 

»-* 

y- 

o 

03 

rn 

OC 

2<X 

It 

o 

y- 

< 

w 

r  * 

• 

UJCM 

vx 

o 

y- 

-J 

It 

o 

00 

•— 

> 

• 

< 

z 

H 

>o 

►—CM 

ocm 

UJ  V 

o 

-J 

UJ 

UI 

UJ 

SG 

o 

wh¬ 

Wh- 

-JO 

N- 

o 

-J 

u 

H. 

IV.. 

om 

n 

0.0 

Q.0C 

o 

03 

o 

t-M 

K 

• 

OCh- 

• 

XUJ 

K 

z 

V 

< 

s 

<o 

►“* 

•-Z 

• 

o 

< 

V 

o 

1-0 

UJ 

00  UJ 

r— 

2 

II 

< 

H 

• 

O 

1- 

o 

-j 

UJ 

o 

W<M 

oc 

1-0 

>- 

*— 

z 

X 

< 

00»-i 

oc 

o 

in 

UJ 

-J 

h- 

o*— 

Zl- 

UJ 

o 

• 

oc 

-J 

G 

< 

It 

* 

ocm 

NlUJ 

II 

u 

CO 

h- 

UJ 

o 

03 

a 

y- 

UJN* 

<z 

> 

oo 

o 

o 

O' 

I- 

z 

o 

‘ 

< 

Z(T 

0*1 

o 

z 

o 

o 

r- 

Jf 

L- 

< 

UJ 

(—1 

<* 

1— 

- J 

UJ  • 

<rx 

o 

< 

2 

o 

o 

• 

< 

* 

CM 

z. 

K 

UJ 

o 

CO 

oc 

»- 

• 

o 

II 

CM 

o 

_J 

H 

c$ 

-J 

o 

o 

zv 

1- 

•— 

• 

00 

00 

o 

*n 

o 

-J 

1- 

»-*- 

Ouj 

o 

-3 

h- 

>- 

o 

i— - 

oc. 

G 

o 

HO 

2 

UJ 

o 

II 

Or 

o 

X 

a 

y- 

u 

z- 

com 

oc 

y- 

i-oc 

V 

X 

oc 

o 

II 

UJ 

< 

iO 

< 

UJ 

o— 

V 

o< 

Cl 

00 

CO 

> 

oc 

oc 

oc. 

1-0 

2 

<11- 

•V 

00 

II 

> 

o 

UJ 

z 

CT' 

UJ 

UI 

< 

o 

o 

o 

X 

bT) 

h: 

z 

oco 

UJ 

c*z 

UJ 

*c 

X 

UJ 

2 

o 

CO 

r* 

n 

XCM 

>ao 

UJ 

o 

t3 

V 

m 

V 

fa 

*— « 

5  • 

V 

oo 

UJX 

*— « 

X. 

I- 

M 

o 

o 

1- 

o 

occr 

02 

1- 

Ul 

V 

o 

2 

o 

o 

>- 

o 

V 

a 

ujct' 

»- 

UJ 

o 

V— 

2 

y- 

03 

o 

a 

< 

UJ 

zo 

VQ 

< 

o 

oc 

< 

I- 

o 

o 

o 

oc 

o 

UJ  • 

VUJ 

-J 

►-4 

< 

-J 

Ul 

O 

O' 

03 

CM 

UJ 

oc 

* 

<CL 

o 

I- 

I- 

-J 

2 

• 

CM 

V 

< 

ZIT 

-JCL 

UJ 

z 

»- 

> 

o 

V 

a 

m 

z 

V 

O' 

i— • 

o 

oc 

< 

X 

CO 

z 

3 

O' 

N* 

00 

* 

uv 

UI 

-J 

<c 

UJ 

II 

I- 

cm 

UJ 

o 

O00 

UL 

oc 

00 

o 

-J 

• 

UJ 

• 

00 

UJ 

z 

09' 

UI 

UJ 

o 

o 

oc 

UJ 

00 

u 

1- 

00  LO 

z 

u. 

<c 

CO 

-J/ 

z 

z 

UJ 

«J 

O' 

M 

o 

luz 

< 

00 

X 

oc 

o 

*— « 

< 

y- 

a 

I- 

a 

GO 

oc 

z 

“3 

a 

z 

0C 

UJ 

X 

II 

1- 

II 

oo 

>m 

OI 

I- 

< 

o 

Hw 

00 

< 

UJ 

LU 

-J 

X 

“3 

< 

</) 

00 

UJ 

oo 

<1- 

oc 

oc 

UJ 

CO 

-J 

h- 

< 

l/> 

o 

mJ 

UJ 

UJ 

>- 

G 

oc>o 

>o 

Ul 

»- 

CL 

u 

< 

o 

o 

1- 

o 

oc 

OC 

G 

in 

G 

o 

-J 

wr- 

< 

X 

-J 

z 

z 

< 

o 

< 

UL 

Cl 

LL 

-J 

a 

-J 

0C 

ZQ0 

zoc 

UJ 

UI 

G 

UJ 

N* 

<1 

o 

UJ 

a. 

O 

o 

< 

UJ  t 

<1 

UJ 

> 

UJ  • 

O  uj 

V 

X 

X 

CO 

X 

oc 

00 

oc 

o 

> 

l-CM 

> 

z 

o 

001- 

z 

UI 

-J 

»- 

00 

oc 

oom 

UI 

UJ 

I- 

1Z 

UI 

o 

o 

o 

♦ 

V 

»— « 

o 

o 

o 

h- 

z 

h- 

K 

y- 

o 

o 

oc*-* 

o 

M 

*1 

t—i 

»- 

o 

** 

< 

•J 

UJ 

-J 

UJ 

UJ 

Ul 

z 

•—« 

Cl 

< 

OC 

1- 

y- 

1- 

y- 

< 

o 

1- 

I- 

a 

3 

3 

o 

X 

•J 

M 

H 

CLX 

UJ 

< 

< 

< 

o 

o 

UI 

< 

< 

X 

00 

z 

l/> 

oc 

-J 

* 

y- 

Ih 

a 

0C 

OC 

oc 

CL 

00 

OC 

oc 

0C 

UJ 

M 

UJ 

< 

M 

G 

< 

oo 

oc 

X 

oc 

1- 

03 

-J 

003 

24 


_0ORN-MAVER  ERODED  POTENTIAL  COPPER  LMfTfTiltE  AND  BULLET 

tB.UkLET  (INCOMING  ENCAD*  ==  llQO..  0  EV  T.IME  STEP  -NULT.tP.MER  =  ."ODC  BULLET  INCIDENT  ON  (.TDD)  .FACE 


o 


UJ 

U 


a 


h-  < 


< 

-j  2: 


-  s 

UJ  <x 


oc 

MJ 

Z 


a 


CO 

o  ^ 

O  CM 

10 

O' 


A 


0 

< 

CO 

it 

m 

H 

O' 

CM 

00 

— K 

CM 

-O 

• 

0 

•  * 

• 

• 

<0 

00 

O' 

CO 

00 

o 

CO 


^  o  m 

CM  CM  -3“ 

•  m 

w  *  *  • 

«  <NJ  CVJ 

N  f—  GO 


o  in 

p-  j- 

CD  «0 

-S'  O' 


O'  m 

it  -J 

p-  o 

•  *  • 


cm  p- 


o  in 
•-  in 
p-  co 
•  *  • 
in  co 
cm  o 


co  p- 

o  o 

fO  it 

f*-  o 
CM  O 


CO  it 
fO  CM 
P-  O' 

0s  O' 
cm  in 


O'  p- 

r-  cm 

00  0 

O  N- 

0  0 

0  O' 

•  *  • 

CM 

*— 

CM 

at 

X 


it  o 
CM  *— 


CM 

U 

X 


CD 


•  *  • 

© 


CO 

U 

X 


o  o 
p-  o' 
o  *- 
•  *  • 
m  o 
m  in 


< 

in 

0 

00 

in 

FO 

in 

P-  O' 

0  0 

p- 

0 

0 

m 

it 

r— 

in 

0  p- 

«—  CP 

CM 

*— 

tn 

m 

CM 

CM 

O'  CM 

r-  O 

UL 

•  * 

• 

•  * 

• 

•  * 

• 

•  *  • 

u 

m 

r— 

00 

O 

CM  O' 

«0  CM 

in 

in 

m 

it 

m  m 

m  m 

> 

UJ 

z 

_  it 

Jt 

00 

00 

0 

fn  -3- 

in  *- 

m  ^ 

O' 

in 

*— 

CD 

O 

in  >0 

P-  O' 

*- 

CM 

CM 

0 

O 

P- 

CM  0 

CO  it 

0 

co  •  * 

• 

•  * 

• 

•  * 

• 

•  *  • 

CC 

11  0 

r- 

lit 

nO 

O 

CM 

0  00 

rn  in 

UJ 

tS2  in 

m 

in 

m 

in 

m 

it  CM 

it  CM 

z 

UJ 

0 

N* 

r— 

»— 

O' 

00 

CM 

r— 

00  in 

m  «- 

t- 

O 

m 

*— 

O' 

in 

CM  00 

m  cm 

UJ 

CP 

*— 

m 

it 

<3 

0  0 

O  CM 

z 

•  * 

• 

•  * 

• 

•  * 

• 

m 

m 

rn 

CM 

*— 

0 

0  O' 

in  o' 

¥ 

nO 

CM 

0 

CM 

>0 

CM 

in  r- 

it  r- 

in 

(d 

• 

«— 

P» 

jS 

00 

fO 

in  in 

CM 

O 

0 

00 

h- 

it 

it 

O 

ro  rn 

|Q 

Z 

CM 

CM 

O' 

O' 

0 

CM 

it  — 

X 

•  * 

• 

•  * 

• 

•  * 

• 

•  *  * 

l-> 

ac 

O 

CM 

CM 

r— 

in 

r— 

m  *— 

0 

in 

it 

r— 

r— 

CM  •— 

►— 

O 

CO 

CM 

P- 

O' 

0 

it 

CM 

|| 

co  CM 

O 

m 

P- 

00 

O 

X 

II 

it 

in 

CM 

m 

***’ 

*  •  * 

• 

•  * 

• 

•  * 

• 

N  r- 

m 

co 

in 

it 

in 

s 


o 

CS  N- 

«!*« 

rT  OJ 


0 

CM 

6fi 

X 


9S 


< 


Fig.  6.  TRIANGULAR  COMPUTER  OUTPUT;  (*)  indicates  impact,  point 


Fig.  7.  GEOMETRICAL  EFFECTS  OF  THE  LATTICE  FACE. 
Target  movement  is  more  restricted  if  initial  motion  is  in  shaded  area 
rather  than  in  crosshatched  area. 
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Fig.  8.  GEOMETRICAL  EFFECTS  on  100  ev  bullet;  (*)  indicates  impact  point 


SKSBEC  WENT**!  fcfc'PWE*  LATTICE  ANO  BULLET 

BULLET  INCOMING  ENERGY  «  8C.0  EV  TIME  STEP  MULTIPLIER  »  .C9C  BULLET  INCIDENT  ON  (1001  FACE 


BO  O' 
CD  ^ 
BO  O* 
•  •  • 
m 


A 


H 

< 

H 

O 


o 

B- 

« 

< 

o 

4 

o 

5 

CO 

CO 

O 

4 

CM 

• » 

• 

•  • 

• 

CM 

b- 

o 

o 

m 

B- 

h- 

CO 

b- 

b- 

x 

tii 

_j 

co 

o 

ct 

ui  a 
o 

«-  >• 
b-  o 
*-  o 


2  2 


I-  z 


o 

at 


Ui 

o 

0C 

< 


>  a. 

ui  U 


£  SUJ 


>*  Z 
O  •"* 
at 

Ui  > 

z  o 

Ui  ot 
Ui 

u  z 

**  Ui 
Ui  u 

z  ~ 

**  t- 

*  ui 
z 

to  1-4 


o 

CO 


o 

-4 

B- 

o 

O 

m 

CO 

O 

m 

B- 

BO 

O' 

B*- 

B- 

B- 

m 

O' 

CM 

•  • 

• 

•  • 

• 

«• 

• 

m 

o 

BO 

GO 

B^ 

CM 

o 

B- 

o 

'O 

in 

'O 

in 

4 

o 

CD 

O 

B- 

o 

•— 

O' 

00 

BO 

O' 

BO 

O' 

B~ 

CM 

o 

BO 

O 

B- 

b- 

CM 

•  • 

• 

•  a 

• 

•  •  • 

•  * 

• 

tn 

O 

-4 

CD 

O' 

CM 

4 

tn 

M> 

tn 

in 

4 

in 

4 

4 

00 

m 

O' 

CM 

B- 

O' 

4 

b- 

CM 

CM 

B- 

4 

in 

O' 

O' 

00 

in 

.4 

CO 

BO 

O 

«  • 

• 

•  » 

« 

•  a 

• 

BO 

K 

CM 

in 

CO 

pm 

BO 

4 

-4 

4 

4 

4 

BO 

4 

BO 

B- 

BO 

pm 

CO 

K 

<> 

B- 

CD 

O' 

r- 

O' 

O' 

O. 

o 

tn 

O 

o 

— >  o 

B^ 

o 

«  •-  ' 

•  * 

• 

• 

•  • 

•  • 

• 

*■— 

4 

u 

CM 

CO 

O' 

m 

BO 

BO 

BO 

BO 

BO 

CM 

CM 

CM 

CM 

/ 


/ 


cm  in 

4  O 
4  *- 

♦  a  * 
o  •“ 
m  cm 


* 

•—  4 

CO 

in 

tn 

CM 

o 

BO 

in 

BO 

co 

in 

CM 

in 

o 

O' 

• 

i/> 

4 

CO 

p-  CM 

B- 

'O 

00 

BO 

o 

1-4 

•  •  • 

•  •  • 

• 

•  • 

• 

•  • 

z 

•—  *— 

O 

O' 

o 

B- 

• 

m-  •— 

•o  o 

BO  BO 
O'  B- 
•  •  • 
h-  in 

CM  CM 


CO  O' 
O'  tn 
O'  ro 

•  *  • 

cm  r- 

CM  p— 


o  o 

B-  O 

B-  o 

4  4 

O  B» 

B-  B— 

4  O' 

CM  O 

o'  b- 

in  •— 

in  bo 

*—  ^ 

•  #  • 

M  • 

•  •  • 

•  •  • 

o*  o 

CD  00 

B-  4 

00  o 

•-  CM 

a  o 
o  z 


u 

4  B- 

CM 

O' 

O 

o 

B- 

-o  b- 

B^  B^- 

O' 

B- 

BO 

CM 

CM 

O' 

p-  •“ 

O 

•  •  • 

•  *  • 

•  •  • 

CO 

in 

in 

in 

4 

m 

4 

in  tn 
-a-  in 

O  CM 

•  *  • 

CM  •— 
CM  •— 

(n 

CO 

in 

• 

CM 


CM 

h- 

4  O' 

CM 

BO 

O'  •— 

pm 

oo 

O  B— 

•  m 

• 

•  •  • 

CM 

p*“ 

CM  p- 

o 

in 

• 

M* 

CD 

co 

4 

O 

«*  • 

• 

B) 

28 


Fig.  9.  OUTPUT  TRIANGLE  for  80  ev  bullet;  (*)  indicates  impact  point 
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Fig.  10.  OUTPUT  TRIANGLE  for  40  ev  bullet;  (*)  indicates  impact  point 
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Fig.  II.  OUTPUT  TRIANGLE  for  25  ev  bullet;  (*)  indicates  impact  point 
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Fig.  12.  OUTPUT  TRIANGLE  for  iOQ  ev  ballet;  (*)  indicates  impact  point 
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Fig.  13.  OUTPUT  TRIANGLE  for  80  ev  ballet;  (*)  indicates  impact  point 
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Fig.  14.  OUTPUT  TRIANGLE  for  40  ev  bullet;  (*)  indicates  impact  point 


BORN-MAYER  ERODED  POTENTIAL  COPPER  LATTICE  AND  BULLET 
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Fig.  15.  OUTPUT  TRIANGLE  for  25  ev  bullet;  (*)  indicates  impact  point 


BULLET  INCOMING  ENERGY  =  100.0  EV  TIME  STEP  MULTIPLIER  =  .080  BULLET  INCIDENT  ON  (100)  FACE 
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Fig.  16.  OUTPUT  TRIANGLE  for  100  ev  bullet;  (+)  indicates  impact  point 


BORN-MAYER  EROOEO  POTENTIAL  COPPER  LATTICE  AND  BULLET 
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Fig.  17.  OUTPUT  TRIANGLE  for  80  ev  ballet;  (*)  indicates  impact  point 
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Fig.  18.  OUTPUT  TRIANGLE  for  40  ev  bullet;  (*)  indicates  impact  point 
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Fig.  19.  OUTPUT  TRIANGLE  for  25  ev  bullet;  (*)  indicates  impact  point 


THOMAS-FERMI-FI RSOV  POTENTIAL,  ERODEO,  COPPER-COPPER 
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Fig.  20.  OUTPUT  TRIANGLES  for  5  kev  bullet;  (♦)  indicates  impact  point 


BULLET  INCOMING  ENERGY  =  10000  EV  TIME  STEP  MULTIPLIER  *  .090  BULLET  INCIDENT  ON  (100)  FACE 
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Fig.  21.  OUTPUT  TRIANGLE  for  10  kev  bullet;  (*)  indicates  impact  point 


THO'AS-FERMI-FIRSOV  POTENTIAL*  ERCCEC,  CCPPER-CCPPER 
BULLET  INCOMING  ENERGY  =  30000. C  EV  TIME  STEP  MULTIPLIER  =  .CSC  BULLET  INCIDENT  ON  (ICC)  FACE 
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Fig.  22.  OUTPUT  TRIANGLE  for  30  kev  bullet;  (*)  indicates  impact  point 


BULLET  INCOMING  ENERGY  =5000.0  EV  TIME  STEP  MULTIPLIER  =  .090  BULLET  INCIDENT  ON  (100)  FACE 
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Fig.  23.  OUTPUT  TRIANGLE  for  5  kev  bullet;  (*)  indicates  impact,  point 


THOMAS-FERMI-FIRSOV  POTENTIAL,  ERODED,  COPPER-COPPER 
BULLET  INCOMING  ENERGY  =  10000  EV  TIME  STEP  MULTIPLIER  =  .090  BULLET  INCIDENT  ON  (100)  FACE 
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Fig.  24.  OUTPUT  TRIANGLE  for  10  kev  bullet;  (*)  indicates  impact  point 


THCMAS-'FERM i-fi rsov  pctential,  erccec,  ccpper-ccpper 

BULLET'  INCCMING  ENERGY  =  20000.0  EV  TIME  STEP  MULTIPLIER  =  .CSC  BULLET  INCIDENT  CN  (ICC)  FACE 
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Fig.  25.  OUTPUT  TRIANGLE  for  30  KEV  bullet;  (.)  indicates  impact  point 


APPENDIX  I 


GENERAL  DISCUSSION  OF  THE  PROGRAM 

The  masses  of  the  lattice  atoms  are  equal  to  that  of  copper  and  are 
an  integral  part  of  the  program.  The  bullet  mass  however,  can  be  varied  by 
changing  the  input  data.  The  potentials  and  forces  are  calculated  by  func¬ 
tion  subroutines  at  the  end  of  the  program.  The  main  program  is  not  dis¬ 
rupted  when  a  different  potential  function  is  substituted. 

A  coordinate  system  was  established  with  an  atom  located  at  the  origin. 
For  a  face  centered  cubic  lattice  in  this  coordinate  system,  the  sum  of  the 
coordinates  of  an  atom  is  always  an  even  number.  Fixed  point  numbers  were 
used  for  lattice  units,  and  the  volume  which  contains  all  the  atoms  (movable 
and  fixed)  was  scanned  in  the  x,  y,  and  z  directions.  If  the  sum  of  the 
coordinates  of  a  point  was  even,  an  atom  was  placed  there,  otherwise  a  space 
was  left.  Truncation  of  quotients  in  fixed  point  arithmetic  allows  a  simple 
test  for  odd  or  even  numbers.  The  fixed  point  coordinates  of  the  atoms  were 
stored  in  memory  as  floating  point  numbers  for  later  calculations.  The 
"core"  atoms  were  assigned  numbers  from  2-64  and  the  fixed  atoms  were  number¬ 
ed  from  65  to  172. 

In  the  calculation  of  the  forces  or  potentials  associated  with  an  atom, 
the  atoms  closer  than  nearest  neighbors  in  an  undisturbed  lattice  must  be 
considered.  A  "nested"  Do  loop  was  used  for  these  calculations.  The  out¬ 
side  Do  loop  was  indexed  from  1  to  64  and  the  inner  Do  loops  initial  index 
was  always  one  greater  than  the  outer  Do  and  ended  at  172.  This  arrange¬ 
ment  avoids  the  calculation  of  forces  or  potentials  twice  for  the  same  pair 
of  atoms.  The  x,  y  and  z  coordinate  separations,  in  that  order,  between 
the  atom  under  consideration  (index  of  outer  loop)  and  any  other  atom  were 
tested.  If  any  of  these  separations  were  found  to  be  greater  than  ROE, 
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the  inner  Do  index  was  advanced  by  one.  If  the  coordinates  passed  these 
three  tests,  the  square  of  the  distance  between  atoms  was  compared  to  the 
square  of  ROE  (R0E2).  The  true  distance  of  any  atom  that  passed  this  test 
was  used  in  the  force  or  potential  function  and  the  result  applied  to  both 
atoms.  (In  potential  calculations,  half  the  potential  was  given  to  each 
atom).  If  the  square  of  the  distance  between  atoms  was  equal  to  or  greater 
than  R0E2,  the  inner  Do  index  was  advanced  by  one.  Most  of  the  atoms  will 
be  eliminated  in  the  first  three  tests.  The  procedure  avoids  time  consum¬ 
ing  square  root  calculations  by  the  computer. 

The  double  iteration  method  used  to  move  atoms  requires  two  calcula¬ 
tions  of  resultant  forces  for  each  Mtime  step".  The  same  section  of  the  program 
was  used  for  both  force  calculations,  but  the  equations  required  to  move  the 
atoms  twice  in  the  time  step  are  not  identical.  The  fixed  point  variable 
INDEX  was  used  as  a  switching  device  to  jump  over  the  motion  equations  that 
locate  the  atoms  at  temporary  positions.  The  jump  is  executed  after  the 
second  calculation  of  forces. 

The  movement  of  atoms  is  continued  until  the  kinetic  energy  of  the 
target  is  a  maximum.  It  is  necessary  to  advance  one  time  step  beyond  the 
maximum  to  establish  the  maximum.  The  particle  locations  and  velocities  at 
the  time  of  the  maximum  are  still  in  memory,  and  these  are  used  in  kinetic 
and  potential  energy  calculations  or  as  output. 

Although  the  two  body  problem  is  in  two  dimensions,  the  impact  para¬ 
meter  is  the  same  as  that  in  the  lattice.  The  two  body  problem  is  solved 
in  the  same  manner  as  the  lattice  except  that  the  interaction  is  stopped 
when  the  atoms  are  separated  by  a  distance  equal  to  ROE  which  means  the 
potential  and  force  between  atoms  is  zero.  The  interaction  is  initially 
started  with  the  separation  of  atoms  slightly  less  than  ROE.  The  AT  for 
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the  two  body  problem  is  much  smaller  than  the  one  used  for  lattice  calcula¬ 
tions.  The  small  Zi  T  consistently  produces  very  much  more  accurate  re¬ 
sults  in  the  two  body  problem  than  in  the  lattice. 

The  input  consists  of  the  bullet  energy  and  mass*  its  x  coordinate* 
and  the  time  step  multiplier.  In  order  to  calculate  data  for  all  36  impact 
points,  six  input  cards  are  required;  each  is  identical  except  for  the  x 
coordinate  of  the  bullet.  The  x  coordinates  assigned  are  3.0,  2.9,  2.8, 

2.7,  2.6,  and  2.5  in  that  order.  After  receiving  the  input  from  one  card, 
the  corresponding  z  coordinate  on  the  line  AB  is  calculated  (see  Fig.  6). 

The  computer  proceeds  through  the  program  for  this  point,  then  subtracts  0.1 
lattice  units  from  the  z  coordinate  and  repeats.  It  procedes  down  a  column 
until  all  calculations  for  impact  points  on  the  line  BC  have  been  made.  The 
computer  pauses  at  this  point  to  await  instructions.  Here  it  may  receive  a 
new  data  card  and  start  another  column.  After  the  calculations  for  the  point 
on  the  apex  of  the  triangle  (x*2.5,  zs3a5)  are  complete,  the  computer  can 
be  directed  to  print  the  output  in  the  triangular  form  shown  by  Fig.  6. 
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APPENDIX  II 


BLOCK  DIAGRAM 
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APPENDIX  IE 


DETAILED  FLOW  CHART 


DIMENSION 

AND 

COMMON 

STATEMENTS 


BLOCK  A 


save  memory  locations  for  subscripted 
variables 


L=  \ 


ROE  *  V2 
ROES  *2.0 


impact  point  number 

lattice  unit  distance  between  nearest  neighbors 
and  the  square  of  that  distance 


EPST-  F oRCtffcfe) 


FORCE 

FUNCTION 


gPSTP.'FbTF(ROF)  < 

PoTentiAl 

function! 

_ 

terms  subtracted  from 
calculated  potential  and 
force  that  "erode"  them 


Un9vj  T 


E\/,  BX^  DTI^  PHK3S, 


starts  here  for  every  "column' 
in  impact  area 


?rer  \.o 

-  ls>&  ~  BX 

5  p*(o  =  ex 
^0)  =  e>€ 

PREfcRc- 100.0 
cycLE  =  o.o 

INDex^  o 

ryO)* 

i*  4 

starts  here  for  every  impact 
point;  assigns  coordinates 
to  bullet 
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BLOCK  B 


velocity  components  of 
bullet 


calculates  constants  used 
in  force  and  movement 
equations 


sets  velocity  of  lattice  atoms 
to  zero 
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DO  \qO  J-  1  ,7 


BLOCK  C 

first  numbers  to  be  assigned  to  lattice  atoms 
2  for  movable  core,  65  for  fixed  atoms 


IT,  JT,  KT  are  X,  Y,  Z  coordinates 
in  the  system  used  for  reference 


KT=o 


oosq  k-  ft— (s?) 


l>o 


each  point  with  integer  values  for 
coordinates  is  tested;  if  the  sum  of 
the  coordinates  is  an  even  number, 
an  atom  is  assigned  to  that  position 


DOS8  Is  1,7  * — [SQ\ 


NT C=  IT  -V  JT  t  KT 


these  decision  processes 
determine  if  an  atom  is  in  the 
"core";  this  is  necessary  in 
order  to  have  core  atoms  with 
consecutive  numbers  2-64 
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BLOCK  C  (  con't) 


^  O 
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BLOCK  D 


sets  temporary  and  permanent 
coordinates  of  immovable  atoms 
equal  to  eachother 


sets  initial  potential  and  kinetic  energy 
of  core  atoms  equal  to  zero 

stores  original  positions  of  bullet  and 
core  atoms 
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BLOCK  E 


t; — => 

70 

DO  "15  X- 
FK(X)-0.0 

FY  (X)^O.O 
F?(l)  =  0-° 

/ 

IS 

COkJTItJOE 

initial  forces  between  core  atoms 
set  to  zero 
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BLOCK  E  (con't) 


<0 


PRY=  RY(J)-RY(I) 


RDRY  =  ABSF(DRY) 


DR?  =  R?(J)  -R?(l) 

l  ~ 

ADR2=ABSF(dK*) 


D1S7- 

t»RX8+  DRyV  DR21- 

<f>lST-Ro 

\  X. 

X<o 

D1ST  -  AJdist 

V 

FORCE2  FORC  (M.ST )-fPST  < 

test  Y  coordinate  separation 
between  atom  I  and  J 


test  Z  coordinate  separation 
between  atom  I  and  J 


test  square  of  distance  between 
atoms;  if  square  passes  test, 
actual  distaice  is  found  and  used 
to  calculate  force 


 ^ 

FoftCE 

> 

FuOCTiotJ 

FOD  s  FOfcC^/DlST 

Fft  =  PRX*  FOD 
FB  =  DRY  * 
fC  =  DR?  *  F°D 


reduce  total  force  into  X,  Y,  and  Z 
components 
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BLOCK  E  (con't) 


round  off  procedures  in 
computer  may  lead  to  very 
small  resultant  forces  in 
situations  where  the  forces 
should  be  zero;  this  procedure 
sets  to  zero  any  resultant  force 
less  than  10~*4  newtons 
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BLOCK  FI 
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BLOCK  FI  (con’t) 


VX(1)  =  FX(l)*DTORMl-t*V/K(l) 


Also  SQmfi  -Por  n  £  corwpo/ie^s 


since  mass  of  bullet  can  be 
varied,  equations  of  motion 
for  the  bullet  must  be 
calculated  separately 


aUo  SQIYIC  -for  U  i 

)*DToc  +  fUX(l) 

:  =2  coord  */iq4cs 

_ } 

/ _ 

lNiDEX=l 

V. _ 

SOTO  10 

(in  block e) 


recycles  to  calculate  forces  at 
temporary  locations 
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BLOCK  F2 


'  62 


BLOCK  G 


63 


BLOCK  H 


DIX  =  RlX(9)-MX(l) 
DlY=WY(«)-MY(i) 

DI2=MZ(8)-Mf(l) 


coordinate  separation  between  target  and  bullet 


D81=D1X%  diy2+  DIE* 


distance  between  bullet  and  target 
is  found  and  potential  energy  calculated 


*i»iWD5i 


±_ 


P0Tl6=?cn-F(D8l)-EpSTP" 


PoTeuniftL 

pUfJcTlOU 


I _ . 

PKH1-  0.5i8  *  I08*[\jlx(l)  +  VlY(l)’‘+VJZ(iy1]*  PMflSS 


PKt(l)*PKEl 
t?ot=  o.o 
TPKe  =  o.o 


kinetic  energy  of  bullet 


DO  bOO  1=  I  >4 


lPisI  +  1 


calculate  potential  energy  of  each  atom 
and  the  total  potential  energy 


DObOO  J«  10,112 
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BLOCK  H  (con’t) 


=  ° 


prevents  interaction  of  buLet  with 
immovable  atoms  covering  bombarded 
face 


f<0 

DRY  =  R1Y(J)-R1Y(I) 

f 

/\DRY-ABSF(DRY) 

<^\DRY 

=0  | 

—  KOEy  ^ 

f<° 

DR£  =  R12(J)~  R1H(I) 

J  . . 

Ami-RBSF(ORt)\ 


test  X  coordinate  separation 
between  atom  I  and  J 


test  Y  coordinate  separation 
between  atom  I  and  J 


test  Z  coordinate  separation 
between  atom  I  and  J 


BLOCK  H  (con’t) 
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BLOCK  J 


BLOCK  K 


assign  velocity  to  bullet 
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BLOCK  LI 


68 


69 


BLOCK  R 


70 


BLOCK  W 


1 

^  OUTPUT  \ 

PaeseMTA-now 

OP  OftTM 

v  see  P-iG.  <0  J 


_ &_ 

EM  D 
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BLOCK  Y 


BLOCK  Z 
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APPENDIX  IV 


DEFINITION  OF  VARIABLES 

Rx(i) 

X  coordinate  of  atom  number  I  at  any  time  (lattice  units) 

R1X(I) 

x  coordinate  of  atom  number  I  at  former  position  (same  as 
position  1  in  procedure  section),  stored  while  average 
forces  are  computed  (lattice  units) 

RY(I),  RZ(I),  RlY(I),  R1Z(I)  y  and  z  coordinates  defined  in  the  same 


manner  as  RX(I)  and  RlX(I) 

FX(I) 

x  component  of  force  on  atom  number  I  at  any  time  (newtons) 

FlX(I) 

x  component  of  force  at  former  position  of  atom  number  I 
(newtons) 

FY(I),  FZ(I),  F1Y(I),  F1Z(I)  y  and  z  components  of  force  defined  in  the 


same  manner  as  FX(I)  and  FlX(I) 

VX(I) 

x  component  of  velocity  at  any  time  of  atom  number  I  (m/sec) 

V1X(I) 

x  component  of  velocity  of  atom  number  I  at  former  position 
(m/sec) 

VY(I),  VZ(I),  VlY(I),  V1Z(I)  y  and  z  velocity  components  defined  in  the 
same  manner  as  VX(I)  and  VlX(I) 


C (I) L) 

matrix  for  storing  data  I  from  impact  point  L 

PPE(I) 

potential  energy  of  atom  number  I  (ev) 

PKE(I) 

kinetic  energy  of  atom  number  I  (ev) 

RXR(I),  RYR(I),  RZR(I)  x,  y,  and  z  coordinates  of  all  atoms  in  initial 
positions  (lattice  units) 

ETFP,  AFAC,  ETFF  constants  used  in  Thomas-Fermi-Firsov  potential 


L 

impact  point  number 

ROE 

nearest  neighbors  separation  (lattice  units) 

ROE  2 

ROE  squared 

EPST 

value  subtracted  from  calculated  force  to  erode  the  force 
(newtons) 

EPSTP 

value  subtracted  from  calculated  potential  to  erode  the 
potential  (ev) 
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FORC 

function  subroutine  that  calculates  force;  argument  is 
the  separation  of  atoms  (newtons) 

POTF 

function  subroutine  that  calculates  potential;  argument  is 
the  separation  of  atoms  (ev) 

EV 

original  kinetic  energy  of  bullet  (ev) 

BX 

original  x  coordinate  of  bullet  (lattice  units) 

DTI 

time  step  multiplier 

PMASS 

mass  of  the  bullet  (amu) 

BZ 

original  z  coordinate  of  bullet  (lattice  units) 

BZS 

z  coordinate  of  bullet  at  end  of  column  in  impact  triangle, 
used  as  stopping  point  (lattice  units) 

PKE8R 

kinetic  energy  of  target  at  the  end  of  previous  time  step 
(ev) 

CYCLE 

number  of  time  steps 

INDEX 

switching  device 

DT 

length  of  basic  time  step  (seconds) 

RM 

mass  of  copper  atom  (kg) 

RM1 

mass  of  bullet  atom  (kg) 

DTOC 

a  constant  used  in  movement  equations;  when  multiplied  by 
the  sum  of  two  velocities,  it  computes  the  average,  multiplies 
it  by  DT  and  gives  the  result  in  lattice  units  (sec  x  lattice 
units/m) 

DTORM 

a  constant  used  in  velocity  equations;  it  combines  the 
quotient  DT/RM  into  one  constant  to  avoid  repeated  divisions 
(sec/kg) 

DT02RM 

same  as  DTORM  except  it  averages  the  force 

DTORM1 

same  as  DTORM  but  used  in  bullet  calculations  only,  since 
the  bullet  mass  may  be  different  from  that  of  the  lattice 
atoms 

DT02RH1 

same  as  DT02RM  but  used  in  bullet  calculations  only  since 
bullet  mass  may  be  different  from  that  of  the  lattice  atoms 

M,N 

atom  numbers 

JT 

y  coordinate  of  an  x-z  plane 
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KT 

z  coordinate  of  an  x-y  plane 

IT 

x  coordinate  of  a  y-z  plane 

NTC 

sum  of  x,y,  and  z  coordinates  of  a  point  in  the  lattice 

DRX,DRY9DRZ  difference  between  x,y9  and  z  coordinates  respectively  of 


any  two  atoms  (lattice  units) 

ADRX,  ADRY, 

ADRZ  absolute  value  of  DRX9  DRY,  and  DRZ 

DIST 

distance  between  two  atoms  or  the  square  of  the  distance 
(lattice  units) 

FORCE 

eroded  force  between  two  atoms  (newtons) 

FOD 

FORCE  divided  by  the  distance  between  two  atoms  (newtons/ 
lattice  unit) 

FA,FB9FC 

x,y,  and  z  components  respectively  of  the  force  between  two 
atoms  (newtons) 

PKE8 

kinetic  energy  of  target  atom  in  the  lattice  (ev) 

DIX,DIY,DIZ  x,y,  and  z  coordinate  separation  between  target  and  bullet 


(lattice  units) 

D81 

distance  between  target  and  bullet  in  lattice  or  the  square 
of  the  distance  (lattice  units) 

POT  18 

potential  between  target  and  bullet  (ev) 

PKE1 

kinetic  energy  of  bullet  at  end  of  time  step  (ev) 

TPOT 

total  potential  energy  of  lattice  (ev) 

TPKE 

total  kinetic  energy  of  lattice  (ev) 

POT 

potential  between  any  two  atoms  (ev) 

C0L3 

potential  of  the  bullet  with  respect  to  the  lattice  minus 
the  target  (ev) 

PE8A 

potential  energy  of  the  target  with  respect  to  the  lattice 
minus  the  bullet  (ev) 

TE 

total  energy  of  the  lattice  (ev) 

PERE 

percent  error  in  total  energy  (percent) 

CONV 

time  step  multiplier  conversion  factor  for  two  body  problem; 
makes  two  body  time  step  multiplier  equal  to  0*01  regardless 
of  the  value  used  for  the  lattice  problem 
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DSQ 

square  of  the  impact  parameter  in  the  two  body  problem 
(lattice  units  squared) 

XIP 

impact  parameter  in  the  two  body  problem  (lattice  units) 

RX2.RY2 

x  and  y  coordinates  of  target  at  any  time  in  two  body 
problem  (lattice  units) 

RX2T, RY2X 

x  and  y  coordinates  of  target  at  former  position*  stored 
while  average  force  is  computed  (lattice  units) 

RX1,RY1,RX1T,  RY1T  x  and  y  coordinates  of  the  bullet  in  the  two  body 
problem  defined  in  the  same  manner  as  RX2S RY2„ RX2T,  and 


RY2T 

VX2.VY2 

x  and  y  components  of  target  velocity  in  two  body  problem 
(m/sec) 

VX2T.VY2T 

x  and  y  components  of  target  velocity  in  two  body  problem 
at  former  position  (m/sec) 

VX1,VY1,VX1T,VY1T  x  and  y  components  of  velocity  of  the  bullet  in  the  two 


body  problem  defined  in  the  same  manner  as  VX2, VY2S VX2T* 
and  VY2T 

FX1.FY1 

x  and  y  components  of  force  on  the  bullet  in  the  two  body 
problem  (newtons) 

FX1T,FY1T 

x  and  y  components  of  force  on  the  bullet  in  the  two  body 
problem  at  its  former  position  (newtons) 

FX2, FY2, FX2T, FY2T  x  and  y  components  of  force  on  the  target  in  the  two 
body  problem  defined  in  the  same  manner  as  FXl,  FYlsFXlTs 


and  FY1T 

DX 

x  coordinate  separation  between  atoms  in  the  two  body  problem 
(lattice  units) 

DY 

y  coordinate  separation  between  atoms  in  the  two  body  problem 
(lattice  units) 

P2KE1 

kinetic  energy  of  the  bullet  in  the  two  body  problem  after 
interaction  is  complete  (ev) 

P2KE2 

kinetic  energy  of  the  target  in  the  two  body  problem  after 
interaction  is  complete  (ev) 

TKE 

total  kinetic  energy  in  the  two  body  problem  (ev) 

T2RA 

recoil  angle  in  the  two  body  problem  (degrees) 

B2SA 

scattering  angle  in  the  two  body  problem  (degrees) 
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SRE8 

SREl 

ATLRA 

TLRA 

ABLSA 

BLSA 

RBMT 

EABL 

ET2,ETL, 


resultant  velocity  of  target  in  lattice  (m/sec) 
resultant  velocity  of  bullet  in  lattice  (m/sec) 
quotient  of  two  velocities;  used  to  find  TLRA 
recoil  angle  in  the  lattice  (degrees) 
quotient  of  two  velocities;  used  to  find  BLSA 
scattering  angle  in  the  lattice  (degrees) 
bullet  to  target  mass  ratio 

energy  absorbed  by  the  lattice;  does  not  include  potential 
or  kinetic  energy  of  target  and  bullet  (ev) 

RET,RKEP,RSA2L,RRA2L,RBFI  ratios  used  to  compare  lattice  interaction 
with  two  body  problem;  the  ratios  are  in  terms  of  variables 
already  defined. 
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APPENDIX  V 


PROGRAM  LISTING 


BLOCK  A 

PROGRAM  TFF 

CIMENSION  RX ( 200 ) ,  RY(20C),  RZI20C) 

BIMtNSICN  FX( 200) , FY( 2CC ) Z ( 20C  ) 

DIMtNSIUN  VX ( 2C0 )  ,  VY(2CC)  ,  VZ(2C0) 

OCIMENSI  ON  R1X12C0),  RlY(2CC),  R1Z(200),  FlX(20C),  F1Y(2CC), 

I  F 1 Z ( 200  )  ,  V 1  X ( 200  )  »  V  1 Y ( 2  00 ) ?  V1Z(2CC) 

DIMENSION  C ( 22 »  36  ) , PPE ( 200 ) , PKE ( 200  ) , RXR ( 200 ) , RYR ( 200 ) , RZR(200) 

COMMON  ETFP, AFAC,ETFF 

AFAC=3. 5906 

ETFF=5.9360E-6 

ETFP=6 .6980E+  3 

L=1 

ROE  =  SQPTF  (2.0 

R0E2  =  2.C 

EPST  =  FORC  (  ROE  ) 

EPSTP  =  POTF(ROE) 

10  FORMAT ( 34X54H  THOM AS-F ER M I -F I R SO V  POTENTIAL,  ERODED,  COPPER-COPPER 

3  READ  INPUT  TAPE  2,4,  EV,BX,DTI,  PMAS S 
4  FORMAT (4E 1 0*3 ) 

BZ  =  BX  +  1.0 
BZS  =  6.0-BX 
5  RX(  1)  =  BX 
RZ (1)  =  B? 

PKE8R  =  -  1  00 . r 
CYCLE  =  0.0 
INDEX  =  0 

RY ( 1 )  =  -SQRTF( 2.0 ) 

12  0 FORMAT ( / , 29H  BULLET  INCOMING  LOCATION  X=  F6.3,  5H  Z=  F6.3, 

I I  OH  ENERGY=  F8.1,  24H  TIME  STEP  MULT  IPL I ER=  F5.3,  /  ) 

15  FORMAT ( 1  Hi ) 

BLOCK  B 

V  x  (  1  )  =  0  .  c 

VY (  1  )  =  SQRTF (  EV  »  63. 54/(  1  . CE3  *  PMASS)  )  *  5.511811  E04 
?  VZ( 1 )  *  0.0 

DT  =  DTI/  VY ( 1 )  *  1.807E-1C 
RM=  105.463911  E-27 
DTOC  =  DT/3.614  E-10 
DTORM  =  DT  /  RM 
DT02RM  =  DTORM/2 .0 
RM 1  =  PMASS  »  1 .6 598  E-27 
DTORM 1  =  DT/RM1 
DT02RM1  =  DTORM 1 /  2.0 
DO  250  I  =  2,  64 
VX (  I  )  =  0.0 
VY {  I  )  =  0.0 
250  CONTINUE 
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eiccK  c 


N  =  65 


JT  =  - 

1 

DO  60 

J 

=  1  ,7 

KT  =  C 

DC  59 

K 

=  1,7 

IT  =  C 

DO  53 

I 

=  1,7 

NTC=I T+Jl +KT 

I  F { NTC 

-NTC/ 2* 2) 

5 

2C 

I  F  {  KT 

— 

1  )  3C, 

21 

21 

IF  (  KT 

-5)  22,  2 

?, 

22 

I  F  (  IT 

- 

1)  30, 

2  3 

23 

IF  (  IT 

- 

5)  24, 

24 

24 

IF ( JT- 

C) 

30,  25 

f 

25 

I F  (  J  T  - 

U) 

40,  40 

? 

50 

RX  {  N ) 

= 

IT 

RY  (  N) 

= 

J  I 

RZ  (  N  ) 

= 

KT 

N  =  N 

+ 

T 

1 

GO  TO 

57 

40 

RX { V  )  = 

i  r 

R  Y  {  M  ) 

= 

JT 

R  Z  (  M  ) 

= 

KT 

M  =  M  + 

1 

57 

I T=  I  T  + 

1 

7,20,57 
,  «■:  1 
3C 

,  25 
,  30 
25 


58 

59 

60 


CONTINUE 
KT  =  KT  -* 
CONTI NU£ 
JT  =  JT  + 
CONTINUE 


BLOCK  0 


DO  65  I  =  65,  172 

PPE {  I  )  =  0.C 
R  1  X  (  I  )  =  RX  (  I  ) 

R  1  V  {  I  )  =  R  Y  {  I  ) 

R 1  Z (  I  )  =  RZ { I  ) 

65  CONTINUE 

CO  66  1=  1,64 

PPE {  I  )  =  0  .  C 
PK E ( I  )  =  C.C 
RXR'I  I  )  =RX  (  I  ) 

RYR ( I  )  =  RY {  I  ) 

66  RZR ( I  )  =  RZ (  I  ) 
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PLGCK  E 


70 

DC 

75  I  = 

F  X  ( 

I ) =  C. C 

F  Y  ( 

I )  =  o.c 

FZ  ( 

I)  =  0.0 

75 

CONTI NUL 

DO  200  I  =  1,Mj 
I P 1  =  I  +  1 
DO  20C  J  =  IP1  ,  172 

IF  (  1-1  )  1  30,  1  20,  120 

13C  I  F  (  J-65)  120,  1  *5,  135 

135  I F ( J- 58 )  2C0,  200,  12C 
120  OR X  =  RX ( J )  -  R X (  I  ) 

ADRX  =  A5SF  (  f'R  X ) 

I  F  (  A  0  P  X  -  ROD  UiC,200,2C0 
140  ORY  =  RY( J )  -  KY(  I  ) 

AORY  =  A.1'  SF  (  OR  Y  ) 

I F ( ADRY  -  RCC  )  1^5,200,200 
14  5  DRZ  =  R Z ( J  )  -  RZ (  I  ) 

ADR  Z  =  ARSF(DRZ) 

I F ( AD R Z  -  ROF )  150,200,200 

150  DIST  =  (I  RX  *  DRX  +  CRY « DRY  +  DRZ  *  DR  Z  ) 
IF  (  DISI  -  P.0E2)  155,  2C0,  200 

155  DIST  =  SZ  RTF  (  DIST) 

FORCE  =  FCRC(  DIST)  -  CPST 

FCO  =  FORCE  /  DIST 

FA  =  DRX  •  FDD 

FR  =  ORY  *  FPD 

FC  =  DRZ  *  FDD 

I F (  J-65)  160,165,165 


1  60 

f  x ( j )  =  rx(j)  + 

F  A 

FY ( J )  =  F  Y( J  )  + 

F  R 

fz ( J )  =  rz(J)  + 

FC 

165 

FX  (  I  )  »  F  X (  I  )  - 

F  A 

F  Y  (  I  )  =  f  Y  (  I  )  - 

F  B 

FZ  C  I  )  =  E  Z (  I  )  - 

FC 

200 

COM T I  NOE 

DO  23  5  I  =  1,64 
I  F  (  A  H  s  F  (  F  X  (  I  )  ) 

-  1  . 

ce-14 ) 

2  1C, 

210, 

2  1  5 

2  1C 

FX (  I  )  =  C.O 

215 

I F (  AESF ( FY ( I  )  ) 

-  1  . 

CE-14 ) 

220, 

220, 

225 

220 

FY (  I  )  =  C.O 

225 

I F (  ABSI(FZ(I) 

)  -1 

.CE-14  ) 

23C , 

230 

,  235 

230 

FZ (  I  )  =  C.O 

235 

CONTINUE 

IF  (INDEX)  240 

,  24C,  420 
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FLOCK  FI 


240 


3  OC 


DO  30  C 

I  =  2,  6*1 

R  1  X  {  I  ) 

=  RX  (  I  ) 

211111 

=  RY( I ) 

*  fcZ  U  ) 

fix  n  ) 

=  FX (  1  ) 

F  1  Y  (  l  ) 

=  F  Y  (  1  ) 

fizii) 

=  FZ  (  I  ) 

VI X  (  1  ) 

=  VX ( 1 ) 

V  1  Y  (  I  ) 

=  VY ( I ) 

V  1  Z  (  I  ) 

=  VZ  (.1  ) 

V  X'  (  I  )  = 

F  X  (  I  )  *  DTCRM  + 

VX  {  I  ) 

VY {  I  )  = 

F  Y  (  I  )  «  D  TORR.  + 

VY  (  I  ) 

VZ {  I  )  = 

FZ (  I)  *DT CRR  + 

VZ  {  I  ) 

R  X  (  I  )  = 

(  VXC  I )  +  V  1  X  {  I)  ) 

*  D TCC 

+ 

R  1  X  (  I  ) 

RY (  I  )  = 

{  V Y (  I  )  t  VI Y (  l  )  ) 

»  DTOC 

+ 

R  1  Y  (  I  ) 

K  Z  (  I  )  = 

(  VZ(  I)  +  V 1 Z (  I )  ) 

*  DTOC 

+ 

R  1  Z  {  I  ) 

R 1  X  (  1  ) 
R1 Y(  1  ) 
R 1 Z(  1  > 
FIX  (  1  ) 
F1Y  {  1  ) 
F1Z  (  1  ) 
VI X(  1) 

V  1  Y  (  1  ) 

V  1  Z  {  1) 
VX(  1  ) 

V  Y  (  1  ) 
VZ  {  1  ) 
rx  ( n 
R  Y  (  1  ) 
RZ  (  1  ) 
INDEX 
GO  TO 


7  C 


*  R  X  {  1  ) 

=  RY< 1 ) 

=  RZ  (  1  ) 

=  F  X  {  1  ) 

=  F  Y  (  1  ) 

=  fz  n ) 

=  VX  (  1  ) 

=  vy  t  n 

=  VZ( 1 ) 
FX ( 1 )  « 

FY( 1)  « 

I-  Z  (  1  )  » 
(  V  X  {  1  ) 
(  V  Y  (  1  ) 
{  VZ {  1  ) 
1 


D  TCRR 1 
D  TORH 1 
0  TCRM 
+  V  1  X  (  1  ) 

+  VI Y (1 ) 

+  VIZ  (  1  )  ) 


+  VX {  1  ) 

+  V  Y  (  1  ) 

+  V  Z  (  1  ) 

)  *  DTOC 
)  *D  TCC 

«  p.tcc 


+ 

+ 


R  1  X  (  1  ) 
R  1  Y  (  1  ) 

R  1  Z  (  1  ) 


CLOCK  F2 


420 

DC  425 

I  =  2,64 

• 

V  X  (  I  )  = 

(  f  x  (  n  + 

F  1  X  (  I  )  ) 

*  0  T  C  2  R  M 

+ 

V  1  X  (  I  ) 

V  Y  {  I  )  = 

(  r  y  ( i )  + 

r  i  y  1 1 ) ) 

*  UTC2  RiM 

+ 

V  1  Y  {  I  ) 

VZ (  I  )  = 

(  FZ (  I  )  + 

F  1  Z  (  I  )  ) 

♦0T02RM 

+ 

V  1  Z  (  I  ) 

RX(  I  )  = 

(  VXC  I  )  + 

V  1  X  (  I  )  ) 

«  DTOC 

+ 

R  1  X  {  I  ) 

RY( I )  = 

(  V  Y  (  I  )  + 

V  1  Y  (  I  )  ) 

*  D  TCC 

+ 

R  1  Y  (  I  ) 

RZ(  I)  = 

(  VZ  {  I  )  f 

Vl Z {  I  )  ) 

*  DTOC 

+ 

R  1  Z  (  I  ) 

425 

CONTINUE 
VX(  1  )  = 

( FX(  1  )  + 

F  1  X  (  1  )  ) 

*DT02RN 1 

+ 

V 1  X  (  1  ) 

V  Y  (  1  )  = 

( FY{  1  )  + 

F 1Y (  1  )  ) 

*  DT CPRM  1 

+ 

V  1  Y  (  1  ) 

VZ  (  1  )  = 

(  FZ (  1  )  + 

F  1Z (  1  )  ) 

«  D  T  0  2  R  M  1 

+ 

VI  Z  (  1  ) 

Rx {  1  )  = 

{  VX(  1  )  + 

V1X ( 1 ) 

)  *  DTOC 

+ 

R  1  X  (  1  ) 

R  Y  {  1  )  = 

(  VY (  1  )  + 

V  1  Y  (  1  ) 

)  *DTCC 

+ 

R 1 Y(  1  ) 

R  Z  (  1  )  = 

{  VZ  (  1  )  *■ 

VIZ  (  1  )  ) 

*  DTOC 

+ 

R  1  Z  (  1  ) 

INDEX  = 
CYCLE  = 

0 

CYCLE  + 

1 .0 

FLOCK  G 

PK  EG  =32 

. 91 63268E 

-c e* (vx(C)*vx(8)  + 

V  Y { 8 ) « VY ( 8 ) 

IF(PKtB) 

460,  460 

,455 

455 

IF ( PKE8R 

-PKE8)  460,465,465 

460 

PKE8R  = 
GO  TO  7 C 

PKE8 

• 

465 

CONTINUE 

V  Z  (  8  )  «  V  Z  (  8  )  ) 
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( 


BLCCK  II 


+  D I Z*D I Z ) 


(1 1  X  =  R  1  X  (  R  )  -  R 1  X  (  1  ) 

0  I  Y  =  R  1  Y  (  8  )  -  H  1  Y  (  1  ) 

L'  I  Z  =  R  1  Z  (  8  )  -  R  1  Z  (  1) 

001  =  ( C I  X *  D I  X  +  UIYeCIY 

IF(CP1  -  ROE?)  423,435,4  ... 

428  D81  =  SC'-’TF  (OF  1  ) 

430  POT  1 8  «  P0THD8U*  -  EPSTP 
on  TO  4  4  c 
4  38  POT  18  =  c.O 

4 4C  PKE1  =  C . 5 1804  1C? L-C8* ( V 1 X( 1 ) *V1 X( 1 ) 
1  VI  Z  (  1  )  )*PM ASS 
PKF(  1  )  =  PKE1 
TROT  =  0.0 
TPKE  =  O.C 
no  60C  1  =  1,64 

IP1  =  I  ♦  1 


CO  60 C  J  =IP],  172 
I F (  I  - 1  )  52  5  ,  52  5  ,  54  0 

525  I F (  J  —  6  5 )  540  ,53C,  5  50 

53C  I F ( J  -  88)  60C  ,  6CC.540 
54 C  OR X  =  P.  1  X  (  J  )  -  Rl  X  (  I  ) 

AORX  =  AiiSF(DRX)  ' 

IF  (  At) RX  -  ROF  )  545,600,600 
545  ORY  =  R1Y(J)  -  R1Y(I) 

AORY  =  AH  SF ( DRY ) 

I  F  (  A D R  Y  -  ROE  )  55 '0 , 6  C  C  ,  6 0 0 
550  DR  Z  =  Rl  /  (  J  )  -  Rl Z  (  I  ) 


+  V  1  Y  (  1  )  *  V  1  Y  (  1) 


+  V  1  Z  (  1  ) 


ft 


AF.  R  Z  =  A  f!  S  F  (  0  R  7  ) 

IFtAORZ  -  ROE)  555,600,600 

555  HIST  =  ( l/RX*DRX  +  DR  Y«  DRY  +  CRZ*DRZ) 

I  F ( D I  ST  -  RCF 2  )  5  56  , oOO , 600 

556  DI ST  =  So  RTF ( D I  ST ) 

560  POT  =  PCiF(DIST)  -  EPSTP 
PP E ( I  )  - P P E (  I  )  +  0 .5*  POT 
PPE(J)  =  PPE(J)  +  0.5  *  POT 
T  POT  =  TPOT  +  POT 
600  CONTINUE 

COL  3  =  PP  E ( 1 )  -0.5  »  POT  18 
P!T8A  =  PPE(  8)  -C.5  *  P0T18 
DO  650  1=2,64 

PKE( I  )=32.9l6i2  68E-C8«(VlX(  I ) « V  1  X (  I ) + V 1 Y ( I ) * V  1  Y ( I) +VlZ(  I  ) * V 1 Z (  I  )  ) 
TPKE  =  TPKE  +  PKF (  I  ) 

650  CONTINUE 

TPKE  =  TPKE  +  PKE  1 

TE  =  TPKl.  +  TPOT 

PERE  =  AGSF( (TE-EV)/LV)  *1.CE2 


BLCCK  J 


„  WRITE  OUIPUT  TAPE  3,10 

t  WRITL  OUTPUT  TAPE  3,  12,  8X,  67.,  EV,  DTI 

WRITE  OUTPUT  TAPE  3,  050 
’ 85C  FORMAT (  /,  19H  LATTICE  LOCATIONS  ,/) 

WRITE  OUTPUT  TAPE  3,860 

860  FORMA  T ( 9CH  ATOM  DX  DY  DZ  VX  V Y 

1  VZ  KE  PE  /) 

SHIFT  =  1.CE-C5 
DO  90C  1=  1,64 
E  N I  2  (  C  )  . 

DX  =  R  1  X (  I  ) -RXR (  I  ) 

+  SSK  (OX)  ,  SLJIL  +  2)  ,  +  SC  M  (  77777B  )  »+FSB( SHIFT  )  ,AJP3(865)  ,  IM7(  1  )  . 

3  65  0Y=  R  1  Y  (  i  )  -  R YR(  I  ) 

+SSK ( DY ) ,SLJ( L  +  2)  ,+ SCM ( 777776 > ,  +  FSB ( SH I F T ) , A J P 3 ( 87 5 ) ,  I  M 2 (  1  )  . 

875  0Z=  R  1 Z (  l  )-  R Z R ( I  ) 

+  SSK( CZ) , SL  J ( L  +2 )  , +  SCM( 777770 ) , +F SR ( SH I F T ) , A JP3 ( 89 5 ) , I N I  2 ( 1 )  . 

885  E N A 2 ( C )  ,  A  J P ( 9CC  ) . 

890  WRITE  OUTPUT  TAPE  3 , 89 5 , I , OX , D Y , OZ  ,  V  I  X ( I ) , V 1 Y ( I  )  ,  V  1  Z (  I  )  ,  P K E (  I  )  , 

1  PPL (  I  ) 

90C  CONTINUE 

895  FORMAT ( I10,3F10.5,3E1C.2,2F10.2) 
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HLCCK  K 


CGN  V  = 
DTOC  = 
UTORK 
I  )  I  ("’  2  R  R 

o  t  o  r  ;*  i 


C.C1/CT1 
C  !  or,  s  COMV 
=  F/TCRK  «•  CCN  V 
=  DTcrRM  *  CONV 

=utcryi  *co,nv 


DT02Rf'l  --  DT02R!'/Il  *  CONV 
RX  2  =  C . 0 
RY2  =  C.C 


rSO  =  (  PZ-3  .0)MBZ-3.0)  +  (BX-3.C  )*  (BX-3.0) 

R Y  1  =  SCRIF(OSO) 

X  I  P  =  R  Y 1 

R  X 1 =  SCRTF(ROr2-OSa  -  1.0E-06 

VX 2  =  0.0 

VY2  =  C.C 

V  Y  1  =  C.C 

FX 1  =  C.C 

FX2=  C.O 

F  Y 1 =  C.C 

p  y  p  —  £  q 

VXl  =  -ScRTFIEV  «  C.C635U  /  PK^SS)  *  S.SlieilECL 
INDEX  =  C 


BLOCK  Li 


1  CCS  DX  =  R X  1  -  R  X 2 
DY  =  RY1-  RY2 
DIST  =  C/*DX  +  DY  <  C Y 
IF  (  Cl ST-R0E2)  1  C 1 C  ,  1C10,  1070 

BLOCK  L2 

1C1C  DIST  =  Si>RTF(DISr  ) 

FORCE  =  FPRC(OIST)  *  ERST 
F X 1  =  DX/  DIST  *  FORCE 
FX2  =  -FX 1 

F Y  1  =  DY/DIST  *  FORCE 
FY2  =  -  T  Y 1 

I F (  INDEX)  1015,  10  15,  1050 

BLOCK  Ml 


1015  RX1T  =  Ra 1 
RY1T  =  RY1 
RX2T  =  R X  2 
RY2T  =  R  Y  2 
FX 1 T  =  FX  1 
FY1 T  =  FY  1 
FX2T  =  F  X 2 
FY2T  =  F  Y  2 
VX 1 T  =  Vx  1 
VY1T  =  VY1 
V X 2  T  =  VX  2 

V  Y?  T  =  V Y_ 

V  X  1  =  F  X  i  * 

V  Y  1  =  F  Y  •  * 

VX2  =  F  X_  * 
vy2  =  ty;  * 

RX 1  =  (VXl 

R  Y  1  =  (  V  Y  1 
RX2  =  (VN2 
RY2  =  (VY2 
INDEX  =  1 
GO  TO  1 CC  5 


D 1 CRM 1  +  VXl 
D  TORN  1  +  V  Y 1 

DTORM  +  VX2 
D 1  CRM  +  VY2 
V X 1 T ) «  DTOC 
VY1T)  *  DTOC 
VX2T  >  *  DTOC 
VY2T)  *OTOC 


+  RX  1  T 
+  RY11 

+  RX2T 
+  RY2T 
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BLOCK  M2 


1  C3C  VX 1 =  ( F  X i  +  F  X 1 T )  *  OTC2RM1  +  VX1T 
VY1  =  ( F  Y 1  +  F Y 1 T  )  *  0T02RM1  +  VY1T 

VX 2=  ( FX2  +  F  X 2 T ) *  LT02RM  +  VX2 T 

VY  2  =  (FY.  +  FY2T)  *  0  TU  2  R  M  +  VY2T 

RX  1  =  (  Vx  1  +  VX  IT  )  •*  C TOC  +  RX1T 
RY1  =  ( VY  1  ♦  V  Y  IT  )  *  CTOC  +  RYU 

RX?  =  ( VX 2  +  VX2TJ  *  CTCC  +  RX2T 

RY2  =  (  V  Y  2  +  VY2T.)  *CTCC  +  RY2  T 
INDEX  =  0 
GO  TO  ICC 5 

BLOCK  N 

107C  P2KE1 =C. ' 1 8C4 1 C2E-C9* ! VX  1  *VX1  +  VYT »VY1 ) «PMASS 
P2KE2  =  32  . 9  16  3263 L-08  *!VX2*  VX 2  +  VY2»  V Y^ ) 

TKE  =  P2KE1  +  P2KL2 
WRITE  OUTPUT  TAPE  3,  15 

WRITE  OUTPUT  TAPE  2,10 

WRITE  OUI PUT  TAPE  3,  12,  BX,  BZ,  EV,  DTI 

WRITE  OUTPUT  TAPE  3,  1075 

1075  FORMAT (////,  LCH  RESULTS  OF  SIMPLE  TWO  BODY  INTERACTION  ,/  ) 
WRITE  OUIPUT  TAPE  3,1080,  P2K E 1 , P2 KL2 , TKE , E V 

1 C800FORMAT (  15H  KIN  ENRG  PROJ=  ,  F12.5  ,  18E  KIN  ENRG  T ARG=  , 

1  F12.5,  I 4H  TOT  ENRG=  ,  F12.5,  12H  ENRG  IN=  ,  F1C.1,  /////) 
WRITE  OUTPUT  I A  PE  3,6  5.A 

654  FORMA  T  ( 3 1 H  RESULTS  OF  LATTICE-  INTERACTION  ,/,  ) 

WRITE  OUTPUT  TAPE  3,  655 

655  FORMAT  (  1  <j  H  TARGET  VALUES  ,//,90H  TIME  STEPS  POT  E  TARGET  W/R  TO 
1BULLET  KE  TARGET  POT  E  TARGET  W/R  TO  LA  T -PUL  LE  T  ) 

WRITE  OUTPUT  1  APE  3,660,  CYCLE,  POTlB  ,PKE8R,  P E A3 
660  FORMAT (FI O.C,  1PE25.6, tl8.6,E25.6  ) 

WRITE  OUTPUT  TAPE  3,665,PKE1,  C0L3,  V 1 X (  1  )  ,  V 1 Y (1  )  , V  1  / (  1  ) 
6650FCRMAT!/, 14H  PULLET  V ALU E S , / / , 1 5H  KIN  ENERGY  =  F7.1  ,  29F  POT 

1  E  W/R  TU  LAT-TARGLT  =  F7.1,25H  VELOCITY  IN  M/SEC  X=  1PE9.2, 

2  4H  Y  =  E9.2,  4H  l~-  E9.2  ) 

WRITE  OUTPUT  TAPE  3,  675 

675  FORMAT  (  /,  16H  LATfICE  VALUES  ) 

WRITE  OUTPUT  TAPE  2,  000 

3G00F0RMAT  (  /,  60H  POT  ENERGY  KIN  ENERGY  TCT  ENERGY 

1  PER  ERROR  ) 

WRITE  OUTPUT  TAPE  2,  810,  TPOT  ,  TPKE,  TE,  PERE 
310  FORMAT!  31  1  5.5  ,  F  10. 5  ) 

WRITE  OUIPUT  TAPE  3,  11 CC 

1100  FORMAT!//,  71H  COMPARISON  VALUES  OF  LATTICE  REACH  GN  AGAINST  SIM.PL 
IE  TWO  P 00 Y  PROBLEM  ,/,  65H  BOTH  INTERACTIONS  STOPPED  WHEN  TARGtT  K 
2INETIC  ENERGY  MAXIMUM  ) 
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T2RA  =  Al ANF( AESF ( VY2/VX2 ) )  *57.296 
B2S  A  =  A  I  ANF ( AHSF  ( VV 1 /VX 1  ) )  *  57.296 
SRE8  =  S-.RTFl  PKESR/32.91  6526BE-08} 

SRF1  =  SCRTFIPKEl  /  ( PMAS S* C .  5 1 804 1 C? E-08  )  ) 
ATLRA  =  ABSF(V1Y(  8)  )/SRE8 
I  F ( ATLRA-1  .0)  1105,1110,1110 

11C5  TLRA  =  ACOS  F  (  ATLRA)  *  57.296 


GO  TO  1116 
1110  TLRA  =  C. C 

1115  AP.LSA  =  ABSF(  VI  Y(  1  )  ) /SRE1 

IF ( ARLSA  -1.0)  1116,1118,1118 

1116  DLSA  =  ACOS F (  ARLSA)  *  57.296 
GC  TO  1119 

1  1  1 F  PLS A  =  C.O 

1119  RBMT  =  P ,y  A  S  S  /  63.54 

LT2  =  (P2KL2/EV)* 1 .0E02 
,  ETL  =  (  PKE ; R/EV ) *  1 . CEC2 

(  RET=ETL/LT2 

RK8P=PKC1 /P2KL1 

EARL  =  T^-PKF 1-PKEFR-PPE (8J-PPE ( 1 ) 
RSA2L=B2SA/8LSA 
RRA2L  =  1 2RA/TLRA 
RB F  I  =  PKE1/EV 


WRITE  OUTPUT  TAPE  3  ,  1  1  20  ,  ET2 , ET L, RET, RKEP , RBMT 
1  1 2C0F0RMA  T  (  / ,  4  3H  PERCENT  KL  TRANSFLREO  .  TWO  BCCY  =  F9.3, 

118IH  LATT  I C.  E  =  F9.3,//,  5cH  RATIO  KL  TRANSFER  IN  LATTICE 

2/  KE  TRANSFER  TWO  BOGY  =  F7.6,//,43H  RATIO  KE  PRCJ  IN  LATTICE  / 

3  KE  PPOJ  TWO  CCOY ' =  F11.5,//,  36h  RATIO  BULLET  MASS  /  TARGET  MAS 

4  S  =  F  6 . 5  ) 

WRITE  OUTPUT  TAPE  3  ,  1  1  30  ,  EAB  L  ,  P-2S  A ,  BL  SA  ,  T?  R  A ,  T  LR  A,  RS  A  2  L  ,  RR  A2  L 
1 1 300FORMAT ( / , 33H  POT  +  KE  ABSORBEO  BY  LATTICE  =  F12.5  ,//, 

1  39H  SCATTERING  ANGLE  TWO  BODY  =  F7.2,  1 9H  LAT 

2TICE  =  F  7 . 2 , /  /  ,  3  5  F  RECOIL  ANGLE  TWO  BODY  =  F7.2, 

3  1 9H  L  AT  T  I  Cl  *  F7.2,//,63H  RATIO  SCATTERING  ANCLE  TWC 


4 POD Y  /  SCATTERING  ANGLE  LATTICE  = 
5LE  TWC  BODY  /  RECOIL  ANGLE  LATTICE 
WRITE  OUTPUT  TAPE  5,1140,  XIP,  RB  F I 
\1 400FORMAT (  /,  21 H  IMPACT  PARAMETER  = 

IT  AT  END  OF  RUN  /  ORIGINAL  ENLRGY  = 
WRITE  OUTPUT  TAPE  3, 1 5 


F9.5,//,55H 
F9.5  ) 

F7.4,  54H 

1PE13.5 


RATIO  RECOIL  ANG 


RATIO  KE  BLLLE 
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BLOCK  R 


C( 1 ,L )  =  BX 
C  (  2  ,  L  )  =  hi  Z 
C{ 3,L  )  =  X  IP 

IF {  92  —  3.0)  1  1  60,1  1 6  C ,  1 17C 

I  1 6 C  C(4,L)  =  C.C 

GO  TO  11  PC 

1170  C(L,L)  =  ACCSH (3  Z-3.0) /XIP)  *  57.296 

II  8 C  C  (  5  ,  L  )  =  t T  2*  1  .  OL-O? 

C ( 6 , L  )  =  82  S  A 
C(  7,L  )  =  ETL*  1  .GE-C2 
C ( 8  , L  )  =  BLSA 

C { 9 , L  )  =  RLT 

C( 1 0, L)  =-  PPE (8  )  +  PKE8R 
C (  1  1  ,  L  )  =  TLRA 
C(  12, L)=PKL  1 
C  (  1  i ,  L  )  =  8  2K  E  1 
C(  1L,D=  pkesr 
C( 15,L)=  P2KE2 
C ( 1 6 , L ) =  PERL 
C( 17, L)= CYCLE 
C( IP, L)=  PPE ( 1 J 
C( 19,L)=  PPEI8) 

C ( 20 , L ) =  T2RA 

C( 2 1 , L ) =R1Y( 1 ) 

L  =  L  +  1 

BLOCK  T 

RZ  =  BZ-C.  1 
I  MBZS-BZ  15,5,1200 
1300  CONTINUE 

I F ( SENSE  SWITCH  1 )  3,  LOCO 

4000  CONTINUE 


RLCCK  W 


LK  AX  =  L- 1 

WRITE  OUTPUT  TAPE  3,10 
WRITE  PU 1  PUT  TAPE  3,1350 
135C  FGR I7  A  T  (  /  /  ,  1  1  OH  X 

■  WWo  j.Trthlx  SANGLE 

1  360  WRITE  C  U  I  PUT  TAPE  3  ,  1  4  00  ,  ( C  (  I  ,  J  ) ,  1=  1.11  ) 

1400  FORMAT  (3F  9 . 3  •  F  10.  2,FT0.4,F10.2,FI0.4,Fl0.2,F10.4,F13.4,F10.2) 


Z  6 

‘  M  (  L  )  /  T  M  (  2  ) 


ANGLE  T  M  (  2  ) 

TE  TAR  RANG l 6  ,  /  ) 


CC  1  SCO 


iv=  1,5 


WP.ITC  GUI  PUT  TAPC 
J  =  1  2 
I  =  14 

WRITE  OUTPUT  TAPC 


3,15 


IS  KINETIC  ENERGY  IN  EV  CF  PULLET  IN  LATT 
IS  KINEIIC  ENERGY  IN  EV  OF  TARGET  IN  LATT 


3,1425 

WRITE  OUTPUT  T  A  PL  3,10 
WRITE  OUTPUT  TAPE  3 , 1 450 , EV , 0 T I 
WRITE  OUTPUT  TAPE  3,1450 
14  50  FORMAT  (  /  /  ,  33X5  3  H  TOP  NO. 

1ICE,//,  32X57H  ROT  TOR  NO. 

2  ICE  ) 

150C  WRITE  OUTPUT  TAPE  3 , 1 0 50 ,  (  C ( J , L  )  ,  L  =  1  ,  11  )  ,  (C(I,L),  L=1,U), 

1 ( C ( J , L ) ,  L= 1 2  »  20) ,  (C ( I , L ) ,  L= 1 2 , 20 ) ,  (C(J,L),  L=2 1,27), 

2 ( C (  I , L ) ,  L  =2 1 ,2  7)  ,  ( C ( J  ,  L )  ,  L  =  28,  52),  (C(I,L),  L=2B,32), 

3 ( C ( J , L ) ,  L  =  33 , 35)  ,  (C(I,L),  L  =  33,35),  C ( J , 36 )  ,C ( 1 , 36 ) 

UO  1530  M  =  1,5 
WRITE  GUI  PUT  TAPE  3,15 
J  =  1  2 
1  =  13 

WRITE  OUTPUT  TAPE 
WRITE  OUTPUT  TAPE 
WRITE  OUTPUT  TAPE 
WRITE  OUTPUT  TAPE  5,1510 
1510  FORMAT!//,  33X53M  TCP  NO 
1 T I C E, // , 20X64 H  BOTTOM  NO 
2R0DY  PROBLEM  ) 

1530  WRITE  OUTPUT  TAPE  3, 1 eSO ,  (C ( J,L ) ,  L=l,ll),  CC (  I  , L )  ,  L  =  1 ,  1  1  )  , 

"  (  C  (  I  ,  L  )  ,  L=  1  2 ,  20  ) ,  ( C  (  J  ,  L  )  ,  L=2  1  ,  27  )  , 


3,1425 
3,  10 

3, 1 430 , E V, CT  I 


IS  KINETIC  ENERGY  IN  EV  CF  BULLET  IN  LAT 
IS  KINETIC  ENERGY  IN  EV  CF  PULLET  IN  TWO 


1  (C( J,L) ,  L  = 1 2 , 20 )  , 
2 ( C ( I , L ) ,  L=2 1 , 27) , 


(  C  (  I ,  L  )  ,  L  =2  8 , 3  2  ) , 


( C ( J , L  ) ,  L  =28 , 32  )  , 

3  ( C  (  J ,  L  )  ,  L  =  3  3 , 5  5 )  ,  (C(I,L),  L  =  33,35),  C(  J  ,  36)  ,C  (  1 , 36  ) 

00  1 55C  M  =  1,5 
WRITE  OUTPUT  TAPE  3,15 
J  =  14 
I  =  15 

WRITE  OUTPUT  TAPE  3,1425 
WRITE  OUTPUT  TAPE  3,  1C 
WRITE  OUTPUT  TAPE  3  ,  1  4  30  ,  E  V ,  CT  I 
'  WRITE  O.UIPUT  TAPE  3,1540 

1540  FORMAT!//,  33XS3H  TCP  NO.  IS  KINETIC  ENERGY  IN  EV  CF  TARGET  IN  LAT 
1 T I C  E , / / ,  27X66H  BOTTOM  NO.  IS  KINETIC  ENERGY  IN  EV  OF  TARGET  IN  TW 

j.  2C  8C0Y  PROBLEM  ) 

f 550  WRITE  CUIPUT  TAPE  3 , 1 850 , ( C ( J , L )  ,  L=1,ll),  !C(I,L),  L=l,11), 

1 ( C ( J , L ) ,  L  = 1 2 , 2 0 )  ,  ( C  !  I , L  )  ,  L  =  1 2 , 20 ) ,  (C(J,L),  L =2 1,27), 

2 ( C ( I , L ) ,  L  =  2 1 , 27)  ,  (C ( J , L  ) ,  L=28,52),  (C(I,L),  L=28,32), 

3 ( C ( J , L ) ,  L  =  3  3 , 3  5 )  ,  (C(I,L),  L=33,3S),  C( J , 36)  ,C ( 1 , 36  ) 

WRITE  CUT  PUT  TAPE  3,15 
J=  1  7 
1  =  16 

WRITE  CUIPUT  TAPE  3,1425 
WRITE  CUTPUT  TAPE  3,  1C 
WRITE  CUIPUT  TAPE  5 , 1 4 30 , EV , G T I 
WRITE  CUIPUT  TAPE  3,1560 

1560  FORMA  f  (  /  /  ,  44X32H  TOP  NO.  IS  NUMBER  OF  TIME  S  TEPS,  /  / ,  46X28F  BCTTC.V 
1  NO.  IS  PERCENT  ERROR  ) 

WRITE  OUTPUT  TAPE  3 , 1 350 , ( C ( J , L ) ,  L=1,ll),  !C(I,L),  L=l,ll), 

1 ( C ( J , L ) ,  L-12,20) ,  ( C ( I , L ) ,  L=12,20),  (C(J,L),  L=21,27), 

2 ( C ( I , L ) ,  L-21,27),  ( C ( J , L ) ,  L=28,32),  (C(I,L),  L=28, 32), 

3 ! C ( J , L ) ,  L=33, 35) ,  (C!I,L),  L=33,35),  C ( J , 36 ) , C ( I , 36 ) 
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un  isfic  y=i,2 

WRITE  OUTPUT  TAPI  3,15 
J  =  8 


I  =  6 


WR I  T'  OUTPUT  T  A  PC  3  ,  1  425 
WRITE-  OUTPUT  TAPI-  3,'lC 
W  Pn  I  T  E  OUTPUT  TAPE  3 , 1 4 *0 , EV, OT I 
.  _  WRITE  OUTPUT  TAPE  3,1570 

157C  FORMA  I  {  /  /  ,  40X39H  TCP  fJCJ .  IS  SCATTERING  ANCLE  IN  L A Tl  I C E , / / ,  3 4 X 52H 
1  BOTTOM  \0.  IS  SCATTERING  ANGLO-  OF  TWO  BODY  PROBLEM  ) 

T58C  WRITE  OUTPUT  TAPE  3 ,  1  3 80  , ( C ( J , L > ,  L  =  1  ,  1  1  )  ,  (C(I,L),  L=l,ll), 

1 ( C ( J  ,  L )  ,  1=12,20),  ( C (  I  »  L  )  ,  L  =  1 2 , 20 ) ,  (C(J,L),  L  =2  1,27), 

2 ( C (  I , L  )  »  L-2 1 »  27)  ,  ( C ( J , L  )  ,  1  =  28,32),  ( C (  I  ,  L )  ,  L=28,32), 

3 ( C ( J  ,  L  )  ,  L  =  3  3 ,  2  5 )  ,  (C(I,L),  L=33,35),  C ( J,36 )  ,C ( I , 36) 

00  16C0  M  = 1 , 2 

WRITE  0L1 PUT  TAPE  3,15 

J=1  1 


I  =  2  C 

writ:  output  tape 

WRI TO  CUI PUT  TAPE 
WRITE  OUTPUT  TAPE 
WRITE  OUTPUT  TAPE 
1590  FORMAT!//,  4  2X  25H 
1  TOM  NC.  IS  RECOIL 
1600  WRITE  OUTPUT  TAPE 

1  (  C  (  J  ,  L  )  ,  L=  1 2  »  20) 

2  ( C  (  I  ,  L  )  ,  1  =  21  ,2  7) 

3  ( C  (  J »  L  )  ,  L  =  3  3 , 3  5 ) 
DO  16  20  M=  1 , 5 
WRI TO  OUIPUT  TAPE 
J=  1  9 


2,1  UP  5 

‘,10 

3, 1430, EV, DTI 
3,1590 

TOP  NC.  IS  RECOI.L  ANGLE  IN  L A T T l C E , / / , 3 6 X4 7H  ROT 
ANGLE  OF  TWO  BODY  PROBLEM  ) 

3, 1850, (C( J,L) ,  L-1,11),  ( C  (  I  ,  L )  ,  L=l,ll), 

,  ( L (  I , L  ) ,  L= 1 2 , 2C ) ,  (C(J,L),  L=21,27), 

,  ( C  (  J ,  L  )  ,  L  =  20, 32  ) ,  (  C  (  I  ,  L )  ,  L=2  8,  32), 

,  ( C ( I , L  )  ,  L  =  3 3 , 3 5 ) ,  C( J,i6) ,C( I, 36) 

3,15 


1  =  18 

WRITE  OUTPUT  TAPE  3,  1425 
WRI TO  CLIPUT  TAPE  3,10 
WRITE  OU  T  PUT  TAPE  3 , 1 430  ,  L V  ,  D T I 
WRITE  OUTPUT  TAPE  3,1610 

1610  FORMAT (// ,32X55H  TCP  NO.  IS  POTENTIAL  ENERGY  IN  E V  OF  TARGET  IN  LA 

1 T 1  ICE,//,  31X58H  BCriGM  MO.  IS  POTENTIAL  ENERGY  IN  EV  OF  BULLET  IN 

2  LATTICE  ) 

1620  WRITE  OUTPUT  TAPE  3 , 1  850 , ( C ( J  ,  L ) ,  L=l,ll),  ( C ( I , L)  .  L= 1 ,  1  1  )  , 

1  (  C  (  J,  L  )  ,  L  =  1  2 , 20)  ,  (  C  (  I ,  L  )  ,  L  =  1  2 , 20  )  ,  (C(J,l),  1  =  21,27), 

2 ( C (  I , L  ) ,  L  =  2 1 ,2  7)  ,  ( C ( J , L ) ,  l=2B,52>,  (0(1,1),  L=28,22), 

3  ( C  (  J  ,  L  )  ,  L  =  33,  25)  ,  (C(I,L),  L  =  2-3,35),  C  (  J  ,  56  )  ,  C  (  1 ,  36  ) 

CO  1640  W=1,2 

WRITE  OUTPUT  TAPE  3,15 
J  =  2  1 
1  =  12 

WR I TE  GUI  PUT  TAPE  2, 1 425 
WR  I  TE  OUTPUT  TAPt  3. ,  1  C 
WRITE  OUTPUT  TAPE  3, 1 4 a0 , E V, DT I 
WRITE  OUTPUT  TAPE  3,1650 

<1630  FORMAT! //,  36X47H  TCP  NO.  IS  PULLET  PENETRATION  IN  LATTICE  UNITS,// 
1.32X56H  BOTTOM  NO.  IS  KINETIC  ENERGY  IN  EV  OF  BULLET  IN  LATTICE  ) 
•1640  WRITE  OUTPUT  TAPE  3  ,  1  8  50  ,  ( C  (  J  ,  L  )  ,  L  =  1,U),  (C(I,L),  L=l,11), 

1 ( C ( J , L ) ,  L=  1 2 »  20)  ,  ( C (  I , L  )  »  L= 1 2 , 20 ) ,  (C(J,L),  L=21,27), 

2 (C (  I , L ) ,  L  =  2 1 , 2  7)  ,  ( C ( J , L  )  ,  L  =  28,32),  (C(I,L),  L=28,22), 

3 ! C ( J  ,  L  )  ,  L  =  33 , 35)  ,  (C(I,L),  L  =  33,  55),  C ( J , 36 ) , C ( 1 , 36 ) 

1425  FORMAT (/////////) 

1430  FORMAT!/,  eX25H  RLLLET  INCOMING  ENERGY  =  ,F8.1,31H  lV  TIME  ST 

1EP  MULTIPLIER  =  ,F5.5,  34H  BULLET  INCIDENT  CN  (100)  FACE  ) 

1 35C  FORMAT!////, 3X1  IF  1C. 3,/, 1 1 1H  *  *  * 

'2  «  ,/,3Xl 1F10. 3, //////, 13X9F10.3t/»  18X82F  *  * 

3*  *  *  *  *  *  **/, 

413X9F10. 3, //////, 23X7F10 .3, /, 23X62H  *  *  * 

5*  *  «  *,/,?3X7F 1C. 3,////// , 33X5F10.3 ,/,38X 

642 H  *  *  *  *  *,/, 33X5T 10. 3, ////// ,4 

73X3F1C. 3, /, 4SX22H  *  *  *,/, 43X 3F 1 0. 3 ,//////, 5 5X F  1 0 . 

8  3,/,  56X211  *,/,53XFlC.3> 

ENO 
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BLOCK  Y 


FUNCTION  POTF(DIST) 

COMMON  El FP.AFAC, ETFF 

POTF=ETFP*  (  1.0+ (f.)  IST+AFAC  )**0.803H  )**  (-3.734)  /  DIST 

RETURN 

END 


BLOCK  Z 


FUNCTION  FORC  (CIST) 

COMMON  61 FP,  AFACi ETFF 
Y ■=  (  C I  ST*  AFAC  )  **C.  8C34 
YP= 1 . C  +  Y 
Z=YP*«(-3.724) 

FOR  C  =  E T F F * Z *  (  1  .0+2.9999*Y/YP)  /(DIST*DI  ST) 

RETURN 

END 

END 


r 
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APPENDIX  VI 


POTENTIALS 


1.  Born-Mayer 


A. 


Potential 


4>0)  =  Ae 


rQ  =  nearest  neighbor  separation  at  zero  pressure  and  absolute  zero 
r^  for  copper  =  2.51A 

The  three  different  sets  of  constants  used  in  the  potential  form  are 

16 

identical  to  those  used  by  Gibson  et.  al. 


Potential 

A(ev) 

1 

0.0392 

16.97 

2 

0.051 

13.00 

3 

0.1004 

10.34 

Potential  No.  2  was  used  extensively  in  the  program.  It  can  be  reduced 


to  a  simplified  form  in  the  following  manner: 


-  13.0  (r- a.SS  \)A,SS\ 

o.os i  a  e 


if  distances  are  expressed  in  lattice  units 

(13.0-^.200  r) 

(J)(r)=  0.051  x  e 

in  the  potential  function  subroutine  this  becomes 

POTF  =  0.051  *  EXPF(13„0-DIST  *  9.208) 
Potentials  one  and  three  become; 

1)  POTF  =  0.0392  *  EXPF(16. 97-DIST  *  8.9996) 
3)  POTF  =  0.1004  *  EXPF(10. 34-DIST  *  7.3115) 
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B.  Force 


Sr 


substitution  of  constants  for  Potential  No«  2 

_  _  (l3.o)(o,OS|)  /G“* 

^  r  ~  2.  SS  l 


lV0(r-?.S5l)/j  S5-, 


if  distances  are  expressed  in  lattice  units  and  electron  volts  changed  to 
joules ; 


(\3.o-  ^.208 r j 

—  0.4lbl4  *  10  x  e 


In  the  force  function  subroutine  this  becomes; 

FORC  =  0.41614E-09  *  EXPF(13.0-9. 208*DIST) 
forces  for  potentials  one  and  three  become; 

1)  FORC  =  0.3127E-09*EXPF(16.97-8.9996*DIST) 

2)  FORC  =  0. 6507E-09*EXPF(10. 34-7. 3115*DIST) 


2.  Thomas-Fermi-Firsov 
A.  Potential 


VM  =  1  E 


K‘ 


TP 


Thomas-Fermi  screening  function 


4W  - 


1  + 


\°<  2 


^3 
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R 

ii 

fti 

£ 

*  O.  8034 

«tj=  3.731 

/  z  V3 

qtp  =  (q^e)  Kqh 

(T  +  <TT«  f 3 

a«=  ^/oe- 

first  Bohr  radius  of  hydrogen 

K-\  = 

-U 

2lS233ev  Qtp=  .8<o0io2LM0  m 

let  x  =  r/a  and  ATFR  =  1/a  (in  lattice  units) 

TF  TF 

then  x  =  ATFR  *  DIST  (DIST  in  lattice  units) 


*«=[' •(*.)“] 


ftTFR  *  DJST] 


let  ATFRjfe^  -  AFAC 

4>(x)  = 


o,go3H 

\  +■  (ft  fac  »  dist; 


_,-3.l34 


V(r) 


Etf  *  (*> ) 

2  x  ATFR  x  DIST 


let  ETFP  =  ETp/  2  *  ATFR 
then 


POT  F  = 


ETFP 

DJST 


\.o  + 


0,0cm : 

(/l  FAC*  diet) 


ai34 


92 


B.  Force 


set  x  =  r/aTp 

^  (V(v)) 

'bsr 

a* 

Vvlwl  

_  <K<)  ,  J_  <j(4>oo) 

^  r  a 

Y *  X  dx 

-  v=(v.r 

then 

d(<i>Cit))  -  -ex',  oc,  <t>(x )  Y_ 

dx  3  X  UY 


^  (  MC*') )  =  ET&  >  drF 
^  r-  *  2rl 


I  +  °<lai3 — 

V  +  I 


<t>00 


if  the  following  substitutions  and  appropriate  changes  in  units  are  made 


where 


F ORC  =  ETFF  «  g 
PI  ST  *  "DIST 


ETFF  =  E__  *  a__  /2 
TF  TF 


1.0  +  P.'WS  *  Y 

yp 


x  =  DIST  *  AFAC  (DIST  in  lattice  units) 

Y  =  (*)°'S02W  YP  =  '-°  +  Y 

.  v- 3,131 

*  =  (YP) 


ETFF  and  AFAC  are  defined  at  the  beginning  of  the  program  -  Y,  YP  and  Z 
are  calculated  in  the  subroutine  each  time  it  is  used. 
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